Spectral function of spinless fermions on a one-dimensional lattice 
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We study the spectral function of interacting one-dimensional fermions for an integrable lattice 
model away from half-filling. The divergent power-law singularity of the spectral function near 
the single-particle or single-hole energy is described by an effective x-ray edge type model. At low 
densities and for momentum near the zone boundary, we find a second divergent singularity at higher 
energies which is associated with a two-particle bound state. We use the Bethe ansatz solution of 
the model to calculate the exact singularity exponents for any momentum and for arbitrary values of 
chemical potential and interaction strength in the critical regime. We relate the singularities of the 
spectral function to the long-time decay of the fermion Green's function and compare our predictions 
with numerical results from the time-dependent density matrix renormalization group (tDMRG). 
Our results show that the tDMRG method is able to provide accurate time decay exponents in 
the cases of power-law decay of the Green's function. Some implications for the line shape of the 
dynamical structure factor away from half-filling are also discussed. In addition, the spectral weight 
of the Luttinger liquid result for the dynamical structure factor of the Heisenberg model at zero 
field is compared with the exact two-spinon contribution. 

PACS numbers: 75.10.Pq, 71.10.Pm 



I. INTRODUCTION 

The single-particle spectral function provides direct in- 
formation about the elementary excitations of a system. 
A recent experiment'^ has demonstrated the possibility 
of measuring the spectral function of strongly interact- 
ing cold atomic gases using momentum-resolved radio- 
frequency spectroscopy.^ This technique can be applied 
to various highly controllable realizations of condensed 
matter systems, motivating theoretical predictions for 
the line shapes which can be measured in such experi- 
ments. 

An example of interaction-induced exotic behavior re- 
flected in the single-particle spectral function is that of 
one-dimensional (ID) fermionic systems. ID metals dis- 
tinguish themselves from higher-dimensional Fermi liq- 
uids by the absence of stable quasiparticles in the low- 
energy regime^^ in which they are generally described 
as Luttinger liquids^ While in a Fermi liquid the main 
effect of weak repulsive interactions on the spectral func- 
tion is to broaden the quasiparticle peak^^ in a Luttinger 
liquid the peak is replaced by an asymmetric power-law 
singularity above the single-particle energy, with an ex- 
ponent which depends on the interaction strengthi^ii The 
case of spin-1/2 fermions, for which two spin-charge sep- 
arated singularities are predicted, is relevant for quasi-lD 
conductors 4 The simpler case of spinless fermions, which 
shall be the subject of this paper, can be realized in fully 
spin-polarized ultracold fermionic gases^ 

The Luttinger liquid result for the spectral function 
is only asymptotically exact in the low energy limit be- 
cause it relies on the approximation of linear dispersion of 
the particle and hole excitations about the Fermi points. 
It is known that this approximation yields the correct 



long distance asymptotic behavior of correlation func- 
tions because band curvature terms in the Hamiltonian 
are formally irrelevant.^ However, recently it has been 
pointed out that dispersion nonlinearity can modify the 
line shape of dynamical correlation functions in the vicin- 
ity of the single-particle energy, u) = e(fc)^iii In fact, 
in the case of weakly interacting spinless fermions with 
parabolic dispersion, the support of the spectral function 
A{k,uj) extends below the single-particle energy. Kho- 
das et alfii have shown that, for generic two-body in- 
teractions, the coupling of the single particle to a con- 
tinuum of excitations with multiple particle-hole pairs 
gives rise to a decay rate and rounds off the singular- 
ity on the single-particle energy. In one dimension, the 
decay rate 1/t is remarkably small because it depends 
on three-body scattering processes. For momentum k 
near and above the Fermi momentum kp, it vanishes 
as 1/t ^ ViiVo - Vk-k^f{k - kpf/m^v^, where 
is the Fourier transform of the interaction potential, vp 
is the Fermi velocity and m is the effective mass. For 
k K. kp, the other energy scale set by band curvature 
is 5uj = {k — kpY /m. For \u} — e(fc)| ~ 1/r, the parti- 
cle spectral function resembles the Lorentzian shape ex- 
pected for a Fermi liquid with quasiparticle decay rate 
1/t. For 1/t <C jw — e(fc)| ^ Sui, the spectral function 
assumes the form of a two-sided power-law singularity 
with a different exponent than the one predicted by Lut- 
tinger liquid theory. The nature of this power law in 
the vicinity of the single-particle energy is understood 
by analogy with the problem of the x-ray edge singular- 
ity in metals.-iS, The power law with the Luttinger liquid 
exponent is only recovered for Slo <^ uj — e(fc) ^ kp/m. 
The crossover between the two power laws is described 
by a universal scaling function of [oj — e{k)]/SujJ^ 



2 



Adding a quadratic term to the dispersion also breaks 
particle-hole symmetry. In general, the hole contribution 
to the spectral function becomes qualitatively different 
from the particle contribution once band curvature ef- 
fects are taken into account. In the continuum model, 
the single hole with momentum —kF<k<kF '\& stable 
since it coincides with the lower threshold of the multi- 
particle continuum. As a result, there is an exact power- 
law singularity at the energy of the single hole. However, 
the exponent in the range |w — e{k)\ <C doj is again not 
the one predicted by Luttinger liquid theoryi^i 

The study of the singularities of spectral functions us- 
ing x-ray edge type effective models is not limited to low 
energies or to weak interactions. In fact, in the case of 
integrable models it is even possible to compute exact ex- 
ponents and energy thresholds for arbitrary momentum 
and interaction strength. The key is to extract these pa- 
rameters from the exact finite size spectrum calculated 
by Bethe ansatz (BA).^'^'^"' The phase shifts that appear 
in the finite size spectrum fix the parameters of the ef- 
fective field theory, which can then be used to calculate 
correlation functions. This was done for the edge singu- 
larities of the dynamical structure factor of the spin-1/2 
XXZ chaiufi^ or equivalently for spinless fermions on a 
lattice,— and for both the dynamical structure factor and 
the spectral function of interacting bosons. The idea of 
extracting exponents of finite-energy spectral functions 
from the BA had appeared earlier in the pseudofermion 
dynamical theory for the ID Hubbard model.— The field 
theory prediction for the singular features of the spec- 
tral function can be combined with numerical methods 
such as the time-dependent density matrix renormaliza- 
tion group (tDMRG)^ to produce high resolution line 
shapes. 

Numerical results for the spectral function of the inte- 
grable model of spinless fermions on a lattice should thus 
offer a quantitative test for the predictions of effective 
x-ray edge models, in particular for the line shape pro- 
posed by Khodas et alfii Two remarks are in order. The 
first one is that integrable models are non-generic in the 
sense that they possess an infinite number of local con- 
served quantities^ What makes these models amenable 
to the BA is precisely that all scattering processes can be 
factorized into a series of two-body collisions. It is not 
clear whether integrable models can have a nonzero de- 
cay rate 1/r which smoothes out the singularities of the 
spectral function. In principle, the question of a small 
versus vanishing decay rate due to interactions can be 
addressed experimentally, since cold atom systems have 
been shown to exhibit very low dissipation and to be ap- 
proximately described by integrable models. 

The second difference from the continuum model is 
that placing the particles in a one-dimensional lattice in- 
troduces effects which are not observed in free space. A 
noteworthy example is the appearance of stable repul- 



sively bound states, which has been demonstrated ex- 
perimentally in the case of bosons3^ Bound molecules of 
spin-polarized fermions have been created artificially by 
sweeping across a p-wave Feshbach resonance.-^ "Anti- 
bound" p-wave molecules in the case of repulsive interac- 
tions should arise naturally as excited states in a lattice. 
An important question is whether these bound states can 
produce a strong response in the single-particle spectral 
function. 

The purpose of this work is to study the spectral 
function for the integrable lattice model of ID spinless 
fermions with repulsive nearest-neighbor interaction. We 
are particularly interested in the regime of high energies 
and strong interactions, where lattice effects are most 
important. 

The paper is organized as follows. In Sec. HIl we dis- 
cuss the kinematics of the model with a cosine dispersion. 
In Sec. IIII[ we present the effective x-ray edge Hamilto- 
nian which describes the behavior of the spectral function 
near the single-particle or single-hole energy. In Sec. IIVI 
we show that repulsively bound states can give rise to 
additional divergent singularities at high energies, which 
are more pronounced at low densities and for momen- 
tum near the zone boundary. In Sec. |Vl we address the 
effects of integrability on the spectral function and on 
the long-time decay of the particle Green's function. We 
show that the broadening of the singularity at the single- 
particle energy in the regime considered by Khodas et al. 
can be recovered by adding an irrelevant interaction term 
to the effective Hamiltonian. In addition, we argue that 
in the integrable lattice model a single particle with low 
velocity can have a nonzero decay rate due to two-body 
scattering processes. In Sec. IVIl we derive the formulas 
for the exact singularity exponents using the BA solution. 
In Sec. IVIIl we compare the analytical predictions for 
the long-time decay of the fermion Green's function with 
high precision data from the tDMRG method. In Sec. 
IVIIIl we discuss the implications of our results for the 
dynamical structure factor away from half-filling. Some 
concluding remarks are offered in Sec. IIXI There are, 
in addition, three appendices. Appendix A contains the 
detailed derivation of the finite size spectrum from BA. 
Appendix B proves the equivalence of our formulas for 
the singularity exponents to those of Ref. JS . Finally, in 
Appendix C we compare the Luttinger liquid result for 
the dynamical structure factor near g = tt to the result 
obtained in the two-spinon approximation, correcting the 
prefactor found in Ref. [l^- 



II. KINEMATICS OF THE LATTICE MODEL 

We consider the model of spinless fermions with 
nearest-neighbor interaction 
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Here, is the operator that annihilates a fermion at site 
j, L is the number of sites, taken to be even, V > is 
the strength of the repulsive interaction, nj = ipji'j is 
the number operator at site j and /z is the chemical po- 
tential. The lowest energy state for A'' fermions is unique 
if we impose periodic (antiperiodic) boundary conditions 
for N odd (even). At half-filling, n = (hj) = 1/2, the 
Hamiltonian is invariant under the particle-hole trans- 
formation ipj — > {—lyi/jj. 

The Hamiltonian of Eq. (1) is equivalent to the 
anisotropic (XXZ) spin-1/2 chain and is BA integrablci^ 
The ground state phase diagram as a function of V and 
fi is known i^°i^^ In this paper we will concern ourselves 
with the region < y < 2 and |^| < 2 -f V, where 
the system is in a gapless phase with power-law decaying 
equal-time correlation functions and low-energy physics 
described by the Luttinger model^. 

The fermion spectral function is defined as 



A{k,Uj) = --lmGret{k,U!), 
TT 



(2.2) 



where 



Gret{k,u) = -i dte'(-+"')*({^,(t),V'i(0)}) (2.3) 



is the retarded single-particle Green's functioui^l Here 
■0^ — e~'^^-'il'j annihilates a fermion with momentum 
k. The spectral function contains both particle and hole 
contributions 



A{k,Lj) = Ap{k,uj) + Ahik,u;), 
which can be written in the form 



(2.4) 



Apik,Lu) = 5{lu~E^+Eo), (2.5) 

a 

Ahik.Lo) = J2\i^\^k\0)^ S{^ + Ea-Eo), (2.6) 



where Eq is the ground state energy and \a) is an exact 
eigenstate of the Hamiltonian (|2.ip with energy Ea. 

In the noninteracting, V = case, the Hamiltonian 
(|2.ip reduces to the tight-binding model: 



(2.7) 



where e(°)(fc) = — 2cos(/c) — /i is the free fermion disper- 
sion and k G (— tt, tt]. The ground state is constructed by 
filling up the single-particle states up to the Fermi level. 



e''°)(fci?) = 0. In terms of fermion density, the Fermi mo- 
mentum reads kp = Trn. Only excited states with a single 
particle with momentum k above the Fermi level couple 
to the ground state via the operator tpl, and contribute 
to Ap{k,uj). As a result, for V = 0: 

4°) (k, uj)^e{\k\~ kF)5{Lo - (fc)). (2.8) 
Likewise, for the hole contribution, 
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f{k,u:)^e{k, 



\k\)6{u 



(fc)). 



(2.9) 



Therefore, in the noninteracting case, the spectral func- 
tion is given by a delta function peak at the energy of 
the single particle (for /sj;- < |fc| < tt) or single hole (for 
\k\ < kp).^ 

In the interacting case, there are nonzero matrix ele- 
ments between the ground state and excited states with 
multiple particle-hole pairs, which then contribute to 
the spectral function. Hereafter we consider the case 
kp < 7r/2 (/i < 0). (From this the case kp > 7r/2 
can be understood by exchanging particles and holes.) 
Since A{k,oj) — A{~k,uj), we also assume < fc < tt. 
Although the renormalized dispersion may deviate from 
the cosine form (unlike the case of parabolic dispersion 
protected by Galilean invariance; c.f. Ref. [Tlh . in this 
section we assume a cosine dispersion in order to discuss 
the support of ^(fc, uj). This is approximately valid in the 
limit of weak interaction. (However, the discussion does 
not depend qualitatively on this assumption; the impor- 
tant feature is simply the existence of a single inflection 
point above the Fermi surface. The exact dispersion can 
be calculated from the BA, as will be done in Sec. I Vill i 

Let us first focus on the particle contribution Ap{k, iv). 
For nonzero interactions, Ap^kjUj) is always nonzero be- 
low the energy of the single-particle excitation. Below 
half-filling, the positive curvature of the dispersion be- 
low the Fermi points implies that the minimum energy for 
fixed total momentum k must correspond to a state with 
|r| -|- 1 particles at the Fermi point with momentum zbfcj?, 
|r| — 1 holes at the Fermi point with momentum ^kp, 

(r) 

and one deep hole with momentum —kp < k\ < kp 
such that 



(2rfcj 



mod 2tt 



(2.10) 



Excitations with a single high energy particle or hole 
can be labeled by three numbers {Nl^Nr, N^)- Nl {Nr) 
is the number of particles created near the left (right) 
Fermi point; Nl,r > denotes particles and Nl,r < 
denotes holes. Nd = 1 in the case of a high energy particle 
and Nd = — 1 in the case of a high energy hole. The 



FIG. 1: Some excitations which define the edges of the sup- 
port of the spectral function: a) excitations with charge +1, 
which contribute to Ap{k,uj); b) excitations with charge — 1, 
which contribute to Ah{k,ijj) (see text for notation). 

excitations that contribute to Ap{k,Lu) have charge +1, 
i.e. Nl + Nii + Nd ~ 1. For instance, the excitations that 
define the deep hole thresholds in Eq. (|2.10[) are labeled 
(-r + l,r+ 1,-1) (see Fig. p. 

Fig. [21 shows the support of the spectral function for 
two rational values of kp/n, kp = 7r/5 and kp/n = 2/5. 
Since in the lattice momentum is only defined mod 27r, 
for every k value there is an infinite number of deep hole 
type excitations, corresponding to different choices of r, 
such that Eq. (|2.10p is verified. Nonetheless, for kp/n 
rational, for general k there is a finite lower threshold for 
the support of Ap{k,uj). The reason is that, for kp/n = 
p/q (Pi 1 integers), the energy of the deep hole thresholds, 

c.(,+i,_,+i,_i)(fc) = -e^'Hkl:^) 

= 2cos ^27r^ - fc^ +^,(2.11) 

is periodic in r with period q. 

If kp/ir is irrational, the energy of the deep hole thresh- 
olds is not periodic in r. This is equivalent to the prob- 
lem of irrational rotations on the unit circle.— In this 
case, for any k there arc infinitely many nondegenerate 
thresholds (— r + 1, r+ 1, —1) and the energy can be made 
arbitrarily low by taking \r\ sufficiently large. In other 
words, for kp/ir irrational, we can move the momentum 
of the deep hole arbitrarily close to the Fermi surface (at 
-iikp) by shifting it by multiples of 2kp (considering only 
the cases where the shift takes the hole to an allowed 
region, below the Fermi surface). Fig. [3] illustrates the 
case kp /-K = !/?> — 5 irrational with (5 <C 1. Consider, for 
example, k = 7r/2. Since kp is not commensurate with 
TT, the energy of the (3, —1,-1) threshold falls slightly 
below the energy of the (0, 2, —1) threshold. Continuing 
with the series, we find that the (6,-4,-1) threshold 
falls at an even lower energy than (3,-1,-1). This se- 
ries of thresholds with decreasing energy implies that the 
support of Ap{k,ijj) extends down to zero energy for all 
k if kp/iT is irrational. 

The edges of the support of the hole contribution 
Ah{k^Lo) can be discussed in a similar fashion. Some ex- 
citations with charge Np^ + Np + Nd = — 1 are illustrated 



,'-'(-3,1,1) '^^.^(-1,-1,1) 



I 1 I 1 I 1 I 1 I 

0.2 0.4 0.6 0.8 




I 1 I 1 I 1 I 1 I 1 I 

0.2 0.4 0.6 0.8 1 

m 

FIG. 2: Support of the spectral function A(k,Lu) for weakly 
interacting fermions with a cosine dispersion (the support is 
the shaded area), for two values of Fermi momentum: a) 
kp = 7r/5; b) kp ~ 2n/5. The positive frequency part corre- 
sponds to the particle contribution Ap{k,ij) and the negative 
frequency part to the hole contribution Ah{k,u>). The lines 
drawn are the dispersion curves for the excitations with multi- 
ple particles and holes at the Fermi points and a single high- 
energy particle or hole, labeled as in Fig. [1] (see also main 
text). The thicker solid line indicates the dispersion of the 
single particle or single hole, lo = ^(k). The other solid lines 
indicate the edges of the support, where A{k, ui) vanishes. 



in Fig. [TJj. The notation is such that the dispersion re- 
lation for the {Np, Npi, Nd) excitation as a function of 
the momentum of the high energy particle or hole is the 
continuation of the dispersion of the {—Npi,—Np,—Nd) 
excitation to negative energies. The upper thresholds of 
Ah(k, Lu) (in the sense of a maximum negative frequencies 
above which AhikjUj) is zero) are defined by deep hole 
excitations labeled (— r, r, — 1). The support oi Ah{k,uj) 
for two cases of rational kp/n is shown in Fig. [2] Note 
that, unlike the case of parabolic dispersion, the single 
hole is not always at the edge of a continuum (see, for 
example, the threshold for kp = 2tt/5 and k < tt/S). 
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FIG. 3: Thresholds for deep hole type excitations for /cf/vt = 
1/3 — 5, 5 ^ 1. The solid line represents the single-particle 
dispersion. For irrational values of ^f/tt and any given fc, 
one can construct a series of "thresholds" (— r + 1, r -|- 1, — 1), 
with increasing number of low-energy particles and holes, 2|r|, 
such that the energy approaches zero. For fc = 7r/2, we have 
'^(0,2,-1) > t^(3,-i,-i) > t^{6,-4,-i) > ■•• • The lower edge 
vanishes for all k and the support of A{k,ijj) covers the entire 
(fc,Cij) plane. 



III. EFFECTIVE HAMILTONIAN FOR 
SINGULARITY NEAR THE SINGLE-PARTICLE 
ENERGY 

The behavior of the spectral function near the energies 
of the single particle or single hole and near the edges 
of the support for nonlinear dispersion can be studied 
using the methods developed in Refs. ^TolfTll. The pro- 
cedure consists of integrating out high-energy modes of 
the fermionic field and introducing subbands which in- 
clude not only the standard low energy modes at ifc^?, 
but also the modes around the momentum of the high 
energy particle (or hole) whose energy defines the thresh- 
old. This leads to an effective model in which the high 
energy particle acts as a mobile impurity coupled to the 



Luttinger liquid modesi^'^'^ The anomalous exponents 
at the thresholds are associated with the phase shifts at 
the Fermi surface due to the creation of the high energy 
particle, in analogy with the x-ray edge singularity in 
metals 

Here we focus on the power-law singularity around the 
single-particle or single-hole energy — e{k). We take 
the continuum limit of the Hamiltonian in Eq. (j2.ip and 
introduce the fermionic field ipi^)- I^i order to study 
the singularity at the energy of the single particle with 
momentum k, we expand 

i^ix) ~ e^'^^^ijRix) + e-'''^^Mx) + e'^^'^x), (3.1) 

where tI'R,l{x) are right and left movers with momentum 
near ±kp and d{x) is the high energy mode defined from 
states with momentum near k, for kp < k < tt. Starting 
from Eq. (|2.ip . we linearize the free fermion dispersion 
e*-°^(/c) about ±kp and k. In the continuum limit, the 
Hamiltonian density for the hopping term in Eq. (|2.ip 
becomes 

Ho = -ivp{ip]^dx^Pi^-^pldxip^)+d\eo~iuodx)d , (3.2) 

where vp — 2 sin fc^^^ is the bare Fermi velocity, Eq — 
e(°)(fc) and uq = 2 sin A: is the bare velocity of the high 
energy particle. Bosonizing right and left movers in the 
form 

i^Rxix) ^ ^_e-'^^«'-(-), (3.3) 

where 4>r,l are the right- and left-moving components of 
a bosonic field and ?7 is a short-distance cutoff, we can 
write 

^0 = Y [(^-'^fl)' + {dM''] +dHeo - tuod^)d (3.4) 

From Eq. (|3.ip . we have the expansion of the density 
operator 



n{x) - n + + i^li^L + + [e'^'""' ^J^J r + e'C^-'^^^^^Vjjd + e'('=+'^'=-)^7/>[d + h.c.]. (3.5) 

Using Eq. (|3.5p and being careful about the point splitting of the fermionic fields, we find the Hamiltonian density 
for the interaction term in Eq. (12. ip 

Tiint ^ — — {dx4>R — dx4>Lf sin fci? (cos A: d^d + isink d'^ dxd ^ 

where we have set 7] — 1 (following Ref. |3i) and omitted irrelevant interaction terms and the renormalization of the 
chemical potential. Note that the coupling to the d modes explicitly breaks the parity symmetry of the original 
Hamiltonian. ^ 

Combining Eqs. p.4p and p.6p . we obtain where (j) and 9 are defined by 

n = ^-^(dxe)\^(dxPj\d^{e~iudx)d 0^^^^^, (3.8) 
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such that [<j){x),dxd{x')\ = i5{x — x' 
V, the parameters in Eq. p.7p are 



To first order in 



V 

VF 

K 



u 



So 



1 



V 



sin kp, 



V 



sin kp, 



(3.9) 

(3.10) 
(3.11) 



The parameters v and K are the famihar renormal- 
ized velocity and Luttinger parameter of the Luttinger 
modeli^ The renormahzed energy and velocity of the 
high-energy particle are given by e and u, respectively. 
Rescaling the bosonic fields by 



(3.12) 



= V/f0 



and introducing the chiral components of the rescaled 
field 



V2 ' 



(3.13) 



we arrive at the effective Hamiltonian density 



n 



where 



1 



+ {dx^hf +S{e-iudx)d 



(klBx^pl - KRdxipFi.)d^d , (3.14) 



kr,l 



l + K 



KRL 



1 -K 



KL,R 



(3.15) 



are the coupling constants between the d particle and the 
Fermi surface modes. 

The Hamiltonian of Eq. (|3.14p is not restricted to weak 
interactions. It can be regarded as exact (up to irrele- 
vant operators) if all the parameters are taken to be the 
renormalizcd ones, to be fixed by a method that is non- 
perturbative in y. In fact, the exact v and K are fixed by 
the low-energy spectrum and compressibility computed 
from the BA solution of Hamiltonian In Sec. Ell 

we will describe how e, u and kr.l can also be fixed with 
the help of the BA solution, so that we end up with a 
parameter-free effective theory for the singularities of the 
spectral function. 

We can decouple the d particle from the bosonic fields 
by performing a unitary transformation 



U = exp 



+ 00 



V27rA' 
with the parameters 

lR.,L{k) 



dx {jRifiR + jL^L)d^d 



(3.16) 



V =F u{k) 

For V 1, we have the weak coupling expressions 
V sin[{kTkF)/2] 



(3.17) 



1R,L ~ T- 



2 cos[(fc± fcF)/2]' 



(3.18) 



In terms of the transformed fields (Pr.l = U^pr^lU^ , d 
UdW, 



dxfRX = dxfR,L ± 



1R,L 

^/2^' 



dU 



d 



d exp 



(3.19) 



(3.20) 



the Hamiltonian (j3.14p becomes nonintcracting (up to 
irrelevant operators) 



^-2 



[dx'-PRf + {dx'-PLf +d'^ {e- iudx) d . (3.21) 



We are now able to calculate the time-dependent particle 
Green's function 



Gp{k,t)^{4,^{t)^lm, 



(3.22) 



using the Hamiltonian of Eq. (|3.2ip . For t ^ e ^, we 
can use the mode expansion of Eq. (j3.ip and write 

Gp{k,t)^ J dx{d{x,t)d^O,0)). (3.23) 

In terms of the decoupled fields, 

Gp{k,t) - / dx{d{x,t)d'' (0,0)) 



X (e 



X e 



-i\/2TTVLVL {x,t) i\/2-KVL^PL (0,0) \ 



,(3.24) 



where 



VR,L{k) 



1 ( i R.Ljk) 
2n 



(3.25) 



are the anomalous exponents. The ground state of 
Hamiltonian (|3.2ip has Nd = high-energy particles, 
and a single d particle is created in the transition. The 
free propagator for the d particle is 



G'(,o)(a;,t) = {d (a;,i)J^(0,0)) 



^5{x-ut). (3.26) 



The correlation functions of free bosonic fields can be 
calculated by standard methods;^ We find 



Gp(fc,i) 



-iet 



[u + v)t — irj 



{u — v)t + irj 

(3.27) 

The long-time asymptotic behavior of the particle 
Green's function is then given by 



Gpik,t) 
where 



exp {—iet — i^l^L ~ sign(u — u)i/fl]} 



v{k) = URik) + VL{k). 



(3.28) 
(3.29) 
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Taking the Fourier transform of Eq. (|3.27p . we obtain 
the particle contribution to the spectral function 



Ap{k,uj) 



+00 



dte"^'Gp{k,t). (3.30) 



The result depends on the sign of u — v. For u > v, 
9{uj — e) sin(7ri'i) + 9{e — uj) sin(7ri/j^) 



Ap{k,uj 
For u < V, 



(3.31) 



A (1 \ ^(^ - e) sin(7r;/) 
Ap{k,uj) — -— . (3.32) 



\uj — e 



Therefore, the effective model predicts that Ap(k,uj) di- 
verges on both sides of the single-particle energy uo — e 
if w > w, but only from above if u < w. In the contin- 
uum model^ the velocity increases monotonically with 
k, thus one always has u > v. In the lattice model, 
the velocity vanishes as k approaches the zone bound- 
aries. In the noninteracting limit, we have u < v for 
TT - kp < k < TT. Note that the result in Eq. p.32p 
does not imply that Ap{k,Lu) vanishes below the single- 
particle energy if u < v. As we discussed in Sec. [TTl the 
support of the spectral function always extends below the 
single-particle energy for kp < 7r/2. However, a singular 
behavior below the single-particle energy only appears if 
the single-particle excitation is unstable (i.e. the energy 
is lowered) against the emission of a low-energy particle- 
hole pair at the Fermi surface pi2 This is the case of a "fast 
particle", with u > v. For a "slow particle" with u < v, 
decay processes with a finite momentum transfer can still 
lower the energy of the single-particle excitation. Both 
finite momentum transfer and three-body collision type 
processes can smooth out the singularity on the single- 
particle energy at the scale of a decay rate 1/t. We shall 
return to this point in Sec. |Vl 

Similar results can be derived for the hole Green's func- 
tion 



G^,{k,t)^{i,l{t)i:,m 



(3.33) 



by employing a Hamiltonian analogous to Eq. (j3.14p and 
interpreting d as the operator that annihilates a hole with 
momentum k < kp. We find 



Gh{k,t) 



exp {-iet - i^[i^L + sign(D - u)i'r]} 



(3.34) 

where £ > and w > are the renormalized energy 
and velocity of the hole, and vr.l are the corresponding 
exponents as given by Eq. (|3.25[) . Typically, we expect 
u < V for a hole below kp < 7t/2, but it is possible that 
the effective mass around the Fermi points changes sign 
in the strongly interacting regime. "^^ For u < v, we find 



Ah{k,-uj) 



{uj - e) s\n[K{vR + vl)] 



(3.35) 



The singularity exponents at the edges of the support 
defined by excitations with multiple particles and holes 
(e.g. (2,0,-1) in Fig. [2]) can also be calculated using 
x-ray edge type effective models as done in Refs. [TTII12I 
In the continuum model, these singularities are always 
convergent due to phase space constraints. In this paper 
we are only interested in divergent singularities, which 
are the most apparent feature of the spectral function 
and govern the long-time behavior of the fermion Green's 
function. Divergent singularities are also the easiest to 
investigate with numerical methods. As we shall discuss 
in Sec. IIVI the lattice model has thresholds defined by ex- 
citations with more than one high energy particle. When 
two of these particles form a bound state, additional di- 
vergent singularities can arise. 



IV. DIVERGENT SINGULARITIES FROM 
BOUND STATES 

The non-monotonic momentum dependence of the ve- 
locity in the lattice model affects the nature of the up- 
per thresholds of multiparticle continua. When an up- 
per threshold is defined by excitations which propagate 
with the same velocity, the interaction between the high- 
energy particles and holes can lead to the formation of 
bound states and strongly affect the line shape of the 
spectral function. In this section, we examine the pos- 
sibility that a two-particle one-hole state can give rise to 
a divergent power-law singularity in the single-particle 
spectral function when two particles with negative mass 
bind for F > 0. 

The existence of bound states above the two-particle 
continuum can be demonstrated directly by solving the 
two-body problem in the lattice. This is a trivial limit 
of the BA equations for the eigenstates of Eq. (I2.1|l for 
N = 2 particles.— Consider the spinless fermion model 
(dropping the chemical potential terms) 

^ = - E(^l^.-+i + '^-^O + Vj2n,n,+i. (4.1) 

3 3 

Consider two-particle eigenstates 



where we may assume 



2J 



The wave function 



must obey the lattice Schrodinger equation 

- - 4>i-i..j - 4>L,3+i - = ^</'ij J (K - i\ > 1) 

4>i+2A - 4>i+i.i-i = - V)(t)i+i^i, (4.3) 

where £ is the energy eigenvalue. We may separate the 
center-of-mass and relative coordinates 



(4.4) 
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Note that the momentum of this state is P. The 
Schrodinger equation becomes 



/(r+l) + /(r-l) 
/(2) 



£ 



2cos(F/2) 

£-V 
' 2 cos(P/2) 



fir) 



/(I)- (4.5) 



The most general solution of the first equation is (for 
r > 0) 



fir) 



(4.6) 



where Q must be either pure real or pure imaginary in 
order for £ to be real. The case of real positive Q corre- 
sponds to a bound state. Q is given by 



cosh Q = — - 



£ 



4cos(P/2)' 

The second of Eqs. (|4.5p determines Q and £ 

V 



= 2 _ 



" 2 cos(F/2) ■ 



(4.7) 



(4.8) 



Thus we see that, for V > 0, a bound state only exists for 
a range of wave vector P where cos(P/2) < 0. Requiring 
Q > 0, we see that 



cos(P/2) < V/2. 



(4.9) 



For small positive V, P/2 = ±{71/2 + 6) with small posi- 
tive 6. Thus the momentum of this state is P « tt ± 2(5. 
The allowed range of momentum is 



sin (5 < V/2. 



(4.10) 



The dispersion relation of this two-particle bound state 



£:(P) =F+| + |cosP, (4.11) 

with the minimum, £ — V, ai momentum P = tt. Note 
that the effective hopping amplitude is tefj — 2/V > 1, 
therefore the repulsively bound state is lighter than the 
free particles. The two-particle continuum in the range 
of Eq. (|4.10p extends over 



4 sin (5 < £ < 4 sin 5. 



(4.12) 



Since 4 sin (5 < -I- (4/y) sin^ 6, the energy of the bound 
state is above the continuum. In this sense, it is an anti- 
bound state (see Fig. 2^). 

The BA solution shows that two-particle bound states 
also exist in the case of finite fermion density.^° We are 
interested in excitations in which a two-particle bound 
state is added to the ground state of the system with 
N particles. The momentum range and dispersion rela- 
tion for this type of excitation can be calculated exactly 
in the thermodynamic limit, as we will discuss in Sec. 
IVIl In the remainder of this section, we work out the 




FIG. 4; (Color online) a) Two-particle continuum and bound 
state band for the two-body problem on the lattice, b) The 
upper threshold of the two-particle continuum is defined by 
two particles with approximately the same momentum and 
negative efi'ective mass. 



effective field theory for the singularities involving two- 
particle bound states. 

We first look more closely at the nature of the upper 
threshold of the two-particle continuum discussed above. 
In general, one could define two subbands with momenta 
around ki and ^2, with ki + k2 = P mod 27r, and expand 
the dispersion relations within each subband in the form 



et{Sh) « £i -f UiSki + 



2m,- 



(4.13) 



for i — 1,2. The energy of a particle pair with total 
momentum P is 



E{5k) = ei{6k) + e2i-5k). 



(4.14) 



Imposing that E{5k) is maximum at 5fc = implies ui = 
U2- That this is a maximum follows from the condition 
TOi — 7712 < 0. This means that, for any dispersion 
relation with a single inflection point above the Fermi 
level, the upper threshold of the two-particle continuum 
for kp < P < has the two particles in the same subband 
with negative effective mass. The wave vectors of the 
particles at the upper threshold are ki = k2 = P/2 — tt 
for kp < P < 7r/2, or ki = k2 = P/2 for 7r/2 < P < tt 
(see Fig. gb). 

Some intuition about the formation of antibound states 
can be gained by considering the case Q ^ 1, so that the 
bound state wave function decays slowly and a continuum 
limit is possible. The Fourier transform of the bound 
state wave function goes as 



Q 



(4.15) 



The bound state wave function involves relative momenta 
of order Q or smaller. Thus it is valid to expand the 
dispersion relation near wave vector k when Q is small. 
We assume we are in this regime and discuss as example 
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the case ki « ^2 
energy field d{x) 



P/2 > n/2. We introduce a high 



d(x) 



(4.16) 



and expand the dispersion in the form e{P/2 +p) « e + 
up + /2m for p ^ 1^*1/2. In contrast with the single- 
particle threshold, we now restrict to the subspace of 
double occupancy of the d subband 



iVd = / dx {d^ (x)d{x)) = 2. 



(4.17) 



It is then important to include in our model the self- 
interaction of the d field. Neglecting for the moment the 
interaction with the Fermi surface modes in the case of fi- 
nite density, we can write down the effective Hamiltonian 
density 



n 



pair 



= dUe- iud^ ~^]d +Vd^dxd d^dU 



2m 



In the weak coupling limit, we have V ~ V and 



1 

771 



V \ P 

1 H smkp cos — - < 0. 
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(4.18) 



(4.19) 



Notice that the interaction between the d particles can- 
not be described by a delta function potential. This is 
because for spinless fermions s-wave scattering vanishes 
identically and the leading process is p-wave scattering. 
In order to obtain a two-particle bound state, we must 
treat the irrelevant interaction term. 

Consider the behavior of the pair spectral function 

dxe~'P^ / dte^'^'npaAx.t), 

J — oo 

(4.20) 

npa»r(a;,i) = (d5,d(a;,i)a,dtdt(o,0)), (4.21) 

near the upper threshold of the two-particle contin- 
uum. For V = 0, the pair spectral function vanishes 
at the upper threshold of the two-particle continuum 
as npair(f, t^) ^ 0{2e - uj)^/2e - lu. For V > 0, low- 
energy scattering between two particles with negative 
mass can give rise to a bound state above this upper 
threshold. The wave function for the two fermionic par- 
ticles can be factorized as ^'(xi, 3:2) ^ e'^^'^"'" f{r), where 
Xcm = {xi+X2)/2 is the coordinate of the center of mass 
and f(r) is the antisymmetric wave function for the rela- 
tive coordinate r = xi — X2- This wave function is found 
by solving the Schrodinger equation 



E-^]f{p) 
m 



J ^Viq)f{p~q), (4.22) 



where V{q) is the Fourier transform of the nearest neigh- 
bor interaction potential (in some approximation with a 
momentum cutoff). 



In one dimension, Eq. (|4.22p with m < always gives 
a normalizable solution with E > for arbitrarily weak 
repulsive interaction. However, the wave function of the 
lowest bound state is an even function of the relative co- 
ordinate. A second bound state with an antisymmetric 
wave function can exist for a potential well with finite 
width and depth if either the interaction is sufficiently 
strong or the particles are sufficiently heavy. According 
to Eq. (|4.19p . the effective mass of the d subband is a 
function of the center-of-mass momentum P and diverges 
as P ^ TT. This corresponds to the momenta of the two 
particles approaching the inflection point of the disper- 
sion (Fig. [4j)). Therefore, for small V > 0, a two-particle 
bound state with a properly antisymmetric wave function 
will form when the two particles have total center-of-mass 
momentum P w tt. Furthermore, the binding energy is 
maximum at P = vr. A bound state dispersion Ei,s{P) 
can be defined for 



iP-vrl < A 



bs 1 



(4.23) 



where hbs ~ OiV) is the half-bandwidth of the bound 
state band. For the integrable model, the exact value of 
h-bs can be obtained from the BA equations (see Sec. |VI] 
B). 

The contribution of a bound state to the pair correla- 
tion function can be written as 



Ilpazr (*^; ^) 



dr 



p 



(4.24) 



In analogy with the definition of d subbands for the 
single-particle dispersion, we take the continuum limit 
in the bound state dispersion and introduce a subband 
with momentum P ^ Pq, with iPo^i"! < ^fcs- We expand 



P-Po+P 



£{P) ~ £bs + UbsP + P^/2mbs, (4.25) 



where Sbs — ^(Po) and Ubs and mbs are the effective 
bound state velocity and mass around momentum Pq. 
Since the bound state dispersion has a minimum at P = 
TT, we have mbs > 0. The propagator for a high energy 
bound state is 



n 



pair 



{x,t) 



dr 



S{x - Ubst). 



(4.26) 



We now turn to the calculation of the contribution 
of the two-bound-particle one-hole state to the single- 
particle Green's function for kp < k < tt. As ex- 
plained in Ref. the Schrieffer- Wolff transformation 
that eliminates from the Hamiltonian off-diagonal scat- 
tering processes (those which take particles out of their 
subbands) generates higher-order contributions to the 
fermionic field. For instance, the leading contribution to 
Gp{k,t) due to a state with two particles in the same 
high-energy subband and one hole at the right Fermi 
point is given by 



Gp{k,t)= J dxe-'''mjix,t)i;''{0,0)), 



(4.27) 
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with 



-ikx 



d\x)d^S{x)^R{x). (4.28) 



The amphtude on the rhs of Eq. is 0{V). The 

idea now is that, if the total wave vector k is in the range 
where the two high-energy particles can bind, the two 
S operators in Eq. (|4.28p can be replaced by a single 
operator corresponding to the creation of a bound state. 
We assume this to be valid both when Q is small and 
when Q is 0(1) (tightly bound pair). 

The minimum momentum fcf, for which a bound state 
contribution to Gp{k, t) exists is obtained by having the 
bound state with the lowest allowed momentum Pq = it— 
h-bs and the hole with the maximum allowed momentum 
kp: 



fch = TT — A 



bs 



(4.29) 



For k > kb, there is a continuum of states defined by a 
single pair and a single hole. For the integrable model, 
the energy of an excited state with total momentum k 
can be written as 



E = £{Po) - e(Po - k), 



(4.30) 



where £{Po) and e(Po — k) are the renormalized energies 
of the bound state and hole, respectively. 

The lower threshold of the one-bound-state one-hole 
continuum is determined by minimizing Eq. (|4.30[) with 
respect to Pq for |Po — t^I < A^^ and |Po~fc| < kp- There 
are three possibilities. 

a) First, for Ubsik + kp) < v, the lower threshold is 
given by a hole at kp and a bound state with momentum 
Pq = k + kp. In this case, we consider the correlation 
function for the operator in Eq. (|4.28p . We define the 
pair creation operator 



BHx)^^J drf{r)d^ 



(4.31) 



The field is bosonic (it creates a p-wave "molecule"). 
We project the fermionic field into the subspace of a sin- 
gle bound state 

iP^x) ^ e~'''''B''{x)'il;R{x). (4.32) 

We can now write down an effective Hamiltonian which 
phenomcnologically describes the coupling of the bound 
state to the low-energy modes 



n2 



BHsbs 



iubsdx)B + ^ {dx^nf + {dxVhf 



{nl'dx^L - ^^tdxVR) . (4.33) 



For the integrable model, all the coupling constants can 
be extracted from the BA as discussed in Sec. IVII B. The 
interaction in Eq. (|4.33p can be diagonalized by the uni- 
tary transformation of Eq. (|3.16p with the appropriate 
parameters 



1r,l 



"R.L 



V =F Ubs 



(4.34) 



The contribution to the particle Green's function be- 
comes 



Gp{k, t) ^ dx 



npajr (3^, t) 



{x — vt + iriYrt (x + vt — irj^L 



(4.35) 



where 



(4.36) 



and Ilpair{x,t) is the bound state part of the pair corre- 
lation function given in Eq. (|4.26[) . As a result, in the 
case Ubs{Po = k + kp) < v, the contribution to Gp{k,t) 
is given by 



Gp{k,t) 



{t-iriY 



(4.37) 



Since 7^^* > for F > 0, one can verify in the weak 



coupling limit that v'j^ -I- < 1 and \ — v'j^- 



0{V). 



As Gp{k, t) decays with an exponent smaller than 1, the 
spectral function acquires a divergent edge singularity 



Ap{k,t) ^ e{oj ~ ebs){oJ - ebs)''''^''- 



(4.38) 



This divergence in the spectral function can be inter- 
preted as an x-ray edge singularity due to the repulsive 
interaction between the bound state and the particles at 
the Fermi surface. 

If V is small, the bound state band is narrow and we 
expect Whs(Po) < v for all tt — Kbs < Pq < it + Kbs- In 
this case, the singularity in Eq. (|4.38p is present for all 

TT — Kbs ~ kp < k < TT + Abs ~ kp. 

b) The second possibility when minimizing Eq. (|4.30[) 
is that, for Ubs (tt -I- Abs) < v and k > it + Abs — kp, it 
is still possible to create one-bound-state one-hole exci- 
tations by creating the hole at a state with momentum 
kh = Pq — k < kp. In this case, the lower threshold of the 
continuum has three high-energy particles. For small val- 
ues of kp — kh, the lower threshold is defined by a bound 
state with the maximum momentum Pq = tt + Abs and a 
deep hole with momentum kh = tt + Abs — k, such that 
the velocity of the deep hole is Uh{kh) > Ubs{TT + Abs). 

As we keep increasing k and decreasing kh, there will 
be a value of hole momentum fc^ at which the hole ve- 
locity Uh{k*f^) equals the maximum bound state veloc- 
ity Ubs(Po = TT + Abs). For k > tt + Abs - k^, min- 
imizing Eq. (|4.30p gives a lower threshold defined by 
a bound state and a deep hole with equal velocities: 
UhiPo -k) = UbsiPo) < V. 

The behavior of the particle Green's function near the 
lower threshold with a bound state and a deep hole can 
be calculated using 



il)^{x) 



'B\x)dl{x), 



(4.39) 



where <i)j(a;) creates a hole. The operators in Eq. (|4.39p 
must be defined after solving the three-body problem in 
which the two particles effectively attract each other (due 
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to the negative mass) and repel the extra hole. We as- 
sume that this has been done and the resulting spectrum 
of the one-bound-state one-free-particle type excitations 
has been parametrized as in Eq. (|4.30p . with no resid- 
ual interactions between the high-energy modes. The 
coupling of these high energy particles to the low-energy 
modes can be described by the effective Hamiltonian den- 
sity 



+dl I Sh - iuhdx 



2mb 

92 



B 



2m J '^^ 



1 



1 



{Ki'd.^L - ^Rd^ifu) B^B 
(43x<PL - n'kd^iPR) d\d^. (4.40) 



The particle Green's function derived from Eqs. (|4.39p 
and is 



Gp{k,t) 



dx- 



(x ~ vt + ir])'^R{x + vt ~ irj^L 



, (4.41) 



where gI^*^' (x, t) is the propagator of the free deep hole 
and v'fii^ ^ 0{V^) arc given by Eq. ((3?25)) with the total 

phase shifts ^r.l = Irl + l due to the bound state 
as well as the deep hole. 

One can verify that, as long as the velocity of the deep 
hole is larger than the velocity of the bound state, Eq. 
(I4.4ip does not lead to a divergence at the lower threshold 
of the one-bound-state one-deep-hole continuum. There- 
fore, for (tt -|- Afts ) < V there is no divergent singularity 
in the wave vector range 11 + Kbs — kp < k < -k + K^s — kf^. 

For k > IT + Kbs — k*f^, we must set Uh = uts in Eq. 
(|4.40p . In this case, replacing both the hole and bound 
state propagators by delta functions with the same ar- 
gument as in Eqs. I|3.26p and (|4.26p would lead to a 
divergence. We must instead calculate these propagators 
by keeping the quadratic term in each dispersion relation. 
This gives 



Upa^r{x,t) ^ y^e™-(--"-*)'/^*,(4.42) 



Gh(x,t) ~ f,-^e^t J^^^mUx-u,Arm_ u 
V 2TTit 



43) 



The integral in Eq. (I4.4ip is then dominated by the re- 
gion X w Ufist. In the limit of large t, we obtain 



Gp{k,t) 



mbsirih 



-i{ets+eh.)t 



2'Kit{mbs + ruh) {t — M/)''n+''i.' 
This leads to a divergence in the spectral function 



(4.44) 



A{k, u) 9{uj - Ebs -eh)iuj-ebs-£h) ^ 



In the weak coupling limit, this is essentially the diver- 
gence of the joint density of states of the two-particle 
bound state and hole with the same velocity. 

In summary, in the case Uf,s(7r -f Af,s) < v, the diver- 
gent edge singularity due to a two-particle bound state is 
absent for tt + Abs — kp < k < n — Abs — k^ but resurges 
for TT — Abs — k'^ < k < TT as SL stronger singularity, with 
exponent close to —1/2 for <C 1. 

c) Finally, the third possibility, when minimizing Eq. 
(|4.30p . is that Ubs{Po) > v for allowed values of Pq < 
TT + Abs- This condition is more likely to be observed 
at lower densities (smaller v) and stronger interactions 
(larger Abs). In this divergent edge singularity 

exists for all k > kb- The nature of the lower thresh- 
old changes at k = Pq — kp, where Pq is such that 
Wbs(-fo*) = "^^ For TT - Abs - kp < k < Pq ~ kp, the 
lower threshold is defined by a hole at fc^ and a bound 
state with momentum k + kp. The corresponding sin- 
gularity is the weaker one given by Eq. (|4.38p . For 
Pq — kp < k < TT, the lower threshold is defined by a 
bound state and a deep hole with the same velocity. The 
edge singularity in this case is described by Eq. (|4.45p . 



V. INTEGRABILITY AND POWER-LAW 
SINGULARITIES 

A. Broadening of the single-particle peak for u > v 



The result in Sec. IIIII predicts that the single-particle 
spectral function exhibits a divergent power-law singu- 
larity at the energy of the single-particle excitation. The 
reason for the exact power-law decay of the particle 
Green's function can be traced back to the ballistic prop- 
agation of the high-energy particle, as expressed by the 
free propagator in Eq. p.26p . If the propagation of the 
d particle is made diffusive, in the sense of a nonzero de- 
cay rate l/r in the particle Green's function, the diver- 
gent singularities in the spectral functions are replaced by 
rounded, Lorentzian-like peaks. In real time, the Green's 
function decays exponentially as Gp{k, t) ^ e^*/'^ at large 

The purpose of this section is to argue that a nonzero 
decay rate l/r in the regime u> v ("fast particle") can 
be recovered within the effective field theory approach by 
adding to Eq. (|3.2ip the following irrelevant interaction 
term which couples the d particle to the bosonic fields 



5Ti = gd'^d dx'PRdxVL, 



(5.1) 



(4.45) 



where g is the coupling constant. In the original fermion 
representation, this corresponds to a six-fermion interac- 
tion term: d^dipl^tpRi^L^L- leads to a three-body scat- 
tering process where the high energy d particle scatters 
simultaneously off two particles slightly below the Fermi 
energy (one near the left branch and one near the right 
branch), producing a final state with the momentum of 
the d particle slightly shifted and two particles produced 
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slightly above the Fermi energy (one near the left branch 
and one near the right branch) . The same interaction was 
considered in Ref. l29|, where it was shown to give rise to 
the finite scattering rate of a mobile impurity coupled to 
a Luttingcr liquid. This type of interaction is allowed by 
symmetry and is generated by bosonization of the inter- 
action term in Eq. (|2.ip to second order 'm V. It is also 
generated by the unitary transformation that eliminates 
the marginal coupling between d and the bosonic fields. 
However, since g is 0{V'^), the correct g cannot be fixed 
in the weak coupling limit by simple bosonization. 

Let us assume that g is present in the effective Hamil- 
tonian and look at the consequences. Consider the effec- 
tive Hamiltonian for a single high-energy particle with 
momentum near k. From Eq. (I3.2ip . the time-ordered 
free propagator for the d particle is 



G(°)(a;,i) = {Td{x,t)d\0,0)) = 9{t)e-''''S{x - ut). 



The Fourier transform is 



G^°\p,uj) = {u} ~ e ~ up + irj) \ 



(5.2) 



(5.3) 



where \p\ < rj^^ <^ k — kp is the momentum within 
the subband. We apply perturbation theory in the g 
interaction. The Dyson equation for the d propagator 
gives the dressed propagator 

Gip,u;) = [u;-e-up^j:{p,Lo)]-\ (5.4) 

where T,(p,Lu) is the self-energy. To second order in the 
coupling constant, we have 



S(p,t^) = -ig^ 



dx ( 



/oo 
dte"^'D^°'>{x,t) 
-OO 



xD^°\x,t)G^"\x,t), 



(5.5) 



where 



1 



27r(ut — irj ^ x)'^ 
The imaginary part of the self-energy reads 



(5.6) 



Im E(p, w) 



dt 

— OO 

i{uj—e—up)t 



[{u — v)t + Z77]^[(w + v)t — ITj] 



,(5.7) 



We can see that, if the velocity of the high-energy par- 
ticle is larger than the Fermi velocity, the integrand in 
Eq. (|5.7p has poles both below and above the real axis 
and Im E is nonzero. This is precisely the condition for 
the existence of the two-sided singularity at the single- 
particle energy. In contrast, if d is a deep hole or a high 
energy particle with u < v, Im E vanishes identically. 
Note also that interaction processes such as d^d (dx'fiR)'^ 
or d^d (dx'fL)'^, in which the d particle scatters only right 



movers or only left movers, do not contribute to the imag- 
inary part of the self-energy. 

Assuming u > v, the decay rate for the high-energy 
particle with momentum near k is 



-Im T,{p,u! — e + up) 



(5.8) 



where rj is the short distance cutoff of the bosonic fields. 
Therefore, this approach yields a cutoff-dependent yet 
nonzero decay rate. It is instructive to consider the scal- 
ing of the decay rate in the low-energy limit k kp. In 
this limit, we have 



V ~ [k — kp)/m. 



(5.9) 



Moreover, the cutoff must scale with the separation be- 
tween the low and high-energy subbands, therefore 



C\k^kp\, 



(5.10) 



where C <C 1 is a constant. Eq. (|5.8|) then simplifies to 



1 Cg^ik^kpf 
T 7rm(2ti)2 



(5.11) 



Finally, in the weak coupling limit g must be regarded as 
the amplitude for a decay process in which the high en- 
ergy particle scatters one right mover and one left mover. 
In terms of the original fermions, the operator in Eq. 
(|5.ip couples a single-particle state to a state with three 
particles and two holes 



(5.12) 



For a generic two-body interaction potential, g ^ 
{a\4>l.\0) must be 0{V'^) and must vanish in the limit 
k kp (no s-wave scattering for spinless fermions). Fur- 
thermore, the coupling constant g has dimensions of 1 /m. 
The combination of parameters that satisfies these con- 
ditions is the one derived in Ref. [ill by perturbation 
theory 



g cx 



VoiVo-Vk^k^) 



(5.13) 



Expanding Vq ~ VQ+V^'q^ /2 for a short-range interaction 
and substituing this form into Eq. ()5.11|) . we find 



1 {VoV^'f{k~kpf 
- ^r-^ 



(5.14) 



Eq. (j5.14p recovers the interaction and momentum de- 
pendence of the decay rate found in Ref. [ll| for small V 
and k kp. In this limit, the decay rate is extremely 
small. However, based on the result of Eq. (|5.8p . we 
expect that the decay rate for V and k — kp of order 1 
can be fairly large in a generic model. This would imply 
that, at high energies (but still in the regime u > v) and 
for strong interaction, the particle spectral function for a 
generic (nonintegrable) ID model exhibits a broad peak 
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near the single-particle energy, not much different from 
the higher-dimensional counterpart. 

However, for an integrable model, it is possible (and 
perhaps probable) that g vanishes identically, since it is 
related to the amplitude of a three-body collision, in the 
sense discussed below Eq. (|5.1I) . Thus we conjecture that 
the irrelevant operator in Eq. (|5.ip is absent in the ef- 
fective Hamiltonian for an integrable model. In fact, the 
exact vanishing of the decay rate to all orders in V and 
k — kp requires the absence of infinitely many more ir- 
relevant interaction terms which couple the d particle to 
the bosonic fields and involve higher powers of dx^PR.L- 
This is only possible by fine tuning of the coupling con- 
stants. That integrability poses nontrivial constraints on 
the coupling constants of some irrelevant operators has 
been demonstrated rigorously in Ref. Here we shall 
justify our conjecture in Sec. IVIII bv examining the nu- 
merical tDMRG results for the particle Green's function 
of the integrable model (|2.ip . 



B. Broadening of the single- particle peak for u < v 

In the "slow particle" regime, where the velocity of 
the high-energy particle is smaller than the Fermi veloc- 
ity, the three-body collision type interaction discussed in 
the previous section yields zero decay rate. This regime 
does not occur in the continuum modelfii The motiva- 
tion to consider three-body collision type processes for 
u > V was that in the continuum model the decay rate 
vanishes at second order in the interaction, making it nec- 
essary to go to 0{V^). In this section we point out that 
in the lattice model the decay rate for a single-particle 
excitation with u < v is nonzero at 0(V'^). This is due 
a two-body collision process in which the slow particle 
scatters off one particle deep below the Fermi surface 
and both final particles are, in general, at high energy 
states above the Fermi surface. Therefore, the broaden- 
ing of the single-particle peak for u < v does not require 
three-body collisions and is present even in the integrable 
model. 

We calculate the decay rate in the limit y <C 1 by 
straightforward perturbation theory in the interaction. 
We write the Hamiltonian in the form H = Hq + Hint, 
where Hq is given by Eq. (|2.7p and 

ki,k2,q 

It follows from Fermi's golden rule that the imaginary 
part of the self-energy to 0{V^) is 

ImS(2)(fc,c.) = -| ^ [y(g)-nfc-p-'?)f 

x^(4"2,)e(-ef)0(e("^,) 
x6{u;-ef\~e^^l^ + e(^y). (5.16) 




FIG. 5: Finite-momentum decay process which gives a finite 
decay rate for the high-energy particle with u < v. 

In contrast with Eq. (|5.12p . the final state \a) con- 
sidered in Eq. I|5.16p has only two particles above the 
Fermi level (in general both of them at high energies) 
and the amplitude (a|V'fc|0) is 0{V). If the dispersion is 
parabolic, Im S*^^) vanishes at the single-particle energy 
due to phase space constraints. However, if the disper- 
sion has an inflection point, it is possible to satisfy both 
momentum and energy conservation such that the decay 
rate is nonzero. In order to calculate 1/r to 0{V'^), we 
can use the noninteracting dispersion e'f'^ = —2 cos k — ^. 
Thus the "on-shell" condition 

4"^+4"^-4"^,-4% = o (5-17) 

is satisfied for 

p = n — k , k + kp — TT < q < k — kp. (5.18) 

This decay process is only kinematically allowed for k > 
TT — kp, the momentum range where u > v in the weak 
coupling limit. It involves a finite- momentum scattering 
process in which the high-energy particle decays as a hole 
is created in the state with momentum tt — fc and another 
high-energy particle is created above the Fermi surface 
(see Fig. [S]). There is a continuum of two-particle one- 
hole excitations, corresponding to different choices of q, 
which are degenerate with the single-particle excitation. 
Importantly, the allowed range of q shrinks to zero in the 
limit kp ^ tt/2 and the decay rate vanishes at half-filling. 

Substituting V{q) = 21^ cos g for the nearest neighbor 
interaction potential in Eq. (|5.16p , we find the decay rate 

- = 21/2 cos2 k r~'' dki , . (5.19) 

T Jf,^ sm fci - sm k 

The explicit result is cumbersome, but the behavior of 
1/t as a function of k is illustrated in Fig. [Hlfor kp — tt/5. 
There is a logarithmic divergence at k — > tt — kp. This 
is an artifact of stopping at second order in perturbation 
theory. The decay process for k — > n — kp involves the 
creation of a hole at the Fermi point with momentum 
kp. We expect that at higher orders in V the density of 
states at the Fermi surface should acquire a power-law 
suppression due to Luttinger liquid physics. As a result, 
1/t should vanish as k approaches the lower threshold. A 
similar discontinuity in the decay rate which is removed 
by treating interactions in the final state has been dis- 
cussed in the context of magnons in the Haldane chain 
in a magnetic field^ 
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FIG. 6: Second-order result for the decay rate of the high- 
energy particle in the wave vector range k > tv — kp, for 
kp = •7r/5. 



The kinematic boundary for this decay process can 
be generalized for finite interaction strength V. Assum- 
ing that the renormahzed dispersion has a single inflec- 
tion point above kp-, the minimum value of k for which 
l/r > is obtained from the condition that Eq. (|5.17p be 
satisfied for the hole at the Fermi point with momentum 
kp 



0. 



(5.20) 



Imposing that Eq. (|5.20p is satisfied for g ^ leads 
to the condition u — v. We conclude that the decay 
rate is nonzero and the spectral function has a broadened 
single-particle peak for momentum k in the regime u < v. 
The precise momentum and interaction dependence of 
the decay rate in this regime is an open question. 

In real time, the corresponding contribution to the par- 
ticle Green's function decays exponentially at large 



Gp{k,t) 



exp{~iet — i'Ki'/2 — t/r) 



(5.21) 



VI. 



EXACT ENERGIES AND EXPONENTS 
FROM THE BETHE ANSATZ 



In order to fix the parameters of the effective model 
in Eq. p.l4p . we need a BA calculation of a physical 
quantity which depends on the coupling constants 
or, equivalently, 7j?,l- First, let us show that the field 
theory approach predicts that the finite size spectrum for 
low-energy excitations in the presence of a high-energy 
particle depends on these coupling constants. Following 
Ref. [stI . we can show that the finite size spectrum for the 
Luttinger model (without the high-energy particle) with 
periodic boundary conditions is 



AE^ + KD^ + n+ 



(6.1) 



where AN is integer, D is integer (half-integer) for N 
even (odd) and n± are integers. AA^ measures the num- 
ber of low energy charge excitations. We have 



AA^ 




(6.2) 



Changing the total number of particles changes the 
boundary conditions of the bosonic field in a finite size 
system)^ Likewise, D is the number of low-energy cur- 
rent excitations, in which particles are transferred be- 
tween the two Fermi points. Hence 



D 



dx (V'ljV'i?, - V'IV'l) 



1 



y/SnK Jo 



dx {d^(pL + dx^PR). (6.3) 



In the presence of a single d particle {{d^d) = 1/L), the 
unitary transformation of Eqs. (|3.16p and (|3.19p takes 



AA^ 
D 



AN 

D - 



_ 7fl + IL 

IR - IL 

AttK 



(6.4) 
(6.5) 



The spectrum of the low-energy particle-hole excitations 
about the Fermi points (n± terms in Eq. (|6.ip ) is not 
modified by the creation of the high-energy particle. We 
conclude that the finite size spectrum in the presence 
of the d particle (including the term of 0(1) from Eq. 
(|3.2ip ) must assume the form 



AE 



2ttv 



1 



AN - 



IR. + 7£ 
27r 



K\D - 



IR - IL 

AttK 



(6.6) 



This is the spectrum of a shifted c = 1 conformal field 
theory jy^^ The parameters Jr^l, which determine the ex- 
ponents in Eq. (|3.25p . can be interpreted as renormalized 
phase shifts at the Fermi points due to the creation of the 
d particle. 

We now proceed to calculate the finite size spectrum in 
the BA, following WoynarovichJ^^ We should first point 
out that there is no correspondence between the fermions 
of the effective field theory and the quasiparticles of BA 
eigenstates. However, we can reasonably expect that the 
power-law singularities of the spectral function, if exist- 
ing, will develop at the energy thresholds of the exact 
spectrum. We can focus on a particular (in practice, 
the simplest possible) class of eigenstates, compute its fi- 
nite size spectrum and compare it with the expression in 
Eq. (|6.6p . This procedure does not require assumptions 
about the form factors of states that contribute to the 
spectral function as written in Eqs. (|2.5p and (12. 6p . 
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A. Single particle excitations 

We start from the Bethe equations for the eigenstates 
of the HamUtonian in Eq. (|2.ip with N < L/2 fermionsi^ 



N 



(6.7) 



fc=i 



The quantum numbers Ij are half-integers for N even 
and integers for N odd. The functions po(A) and 9{X) 
are, respectively, the quasimomentum and two-particle 
scattering phase shift as a function of rapidity A 



po(A) = i\og 
0{X) = t\og 



sinh «/2 + A) 
sinh «/2- A) 
sinh {i( + A) 
sinh {i( — A) 



(6.8) 
(6.9) 



where C = arccos(F/2). A given eigenstate is character- 
ized by a set of rapidities {A^} and the corresponding 
energy is 



N 



where 



eo (A) = 



2 sin^ C 
cosh(2A) — cosC 



(6.11) 



is the bare energy of a particle with rapidity A. 

Consider a state defined by taking {Ij} to be the set of 
all integers (or half-integers) between and /_ , which 
are defined by 



1+ = max {Ij } 



I- = min{/j}- 



(6.12) 
(6.13) 



This is a state with no low-energy particle-hole pairs 
(fij. = 0). The total number of particles is 



N 



IN 



h + l 



(6.14) 



The current carried by an eigenstate is given by the dif- 
ference between the number of right and left movers. It 
can be written as 



2D = h+ I, 



N 



(6.15) 



The ground state for N particles corresponds to the 
choice /+ = — /_ = N/2, with all rapidities Xj in Eq. 
(|6.7p real. To create a single-hole excitation we remove a 
quantum number \Ifi\ < N/2, from the set of /j's. A 
single particle excitation is created by adding one quan- 
tum number Ip with |/p| > A^/2. We distinguish between 
single-particle excitations with real and complex rapidi- 
ties. It follows from Eq. (|6.7p that, for real rapidities 



|Aj| < 00, the quantum numbers Ij are restricted to the 
interval N/2 <\Ij \ < I^o, where^ 



N TT-C fL 



-N 



(6.16) 



As we shall see, this implies a maximum momentum for 
single particle excitations with real rapidities. Note, how- 
ever, that we also get real values oipo{X) and 9{X) in Eqs. 
(|6.8p and (|6.9p (therefore a scattering state of unbound 
quasiparticles) by taking complex rapidities of the form 



A„ 



Re(Ap) 



ITT 



(6.17) 



This is called a negative-parity one-string. This type 
of solution has quantum numbers in the interval I^o < 
\Ip\ < L/2, where 



i =^ 

loo 2 



In the limit of large L, Eq. 
for the density of rapidities p[X)^ 



(6.18) 
becomes an equation 



(6.10) p{X,)^^ 



(A,) + iEA-(A,-Afe) + 



L 



(6.19) 

Here the sum is over all the quantum numbers between 
/_ and /+. The kernel K{X) is given by 



KW = 



d0{X) 
dX 



2sin(2C) 



1 - 2 cos2 C + cosh(2A) ' 



(6.20) 



The last term on the rhs of Eq. (|6.19p is the scattering 
phase shift between the particle at Xj and the high-energy 
particle or hole. In the case of a hole with real rapidity 
Aft,, we have 

<i>h{X) = -KiX-Xh). (6.21) 

In the case of a particle (Ap real or complex), we have 

$p(A) =if(A-Ap). (6.22) 

In the thermodynamic limit, it is convenient to intro- 
duce the shorthand notation for the integral operator K 



K-f{X) 



dx'K{x-x')f{y), 



(6.23) 



where B — A(/+ = N/2) is the Fermi boundary. The 
density of rapidities in the ground state (in the absence 
of the impurity) is given by the Lieb equation 



PGsiX) 



P'oW 
2n 



(6.24) 



The Lieb equation has to be solved consistently with the 
condition for the average density 



dXpGs{X) = n. 



(6.25) 
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The derivation of the finite size spectrum is standard 
and is detailed in Appendix |^ The result for a single 
high-energy particle is 



AS = e(Ap) + — 



(AA^-»Lp)^ 
4K 



K D 



'1771 p ) 



(6.26) 

Here e(A) is the dressed energy defined by the integral 



equationii 



1- J) •e(A)=.o(A), 



(6.27) 



and V is the renornialized velocity given by 

27rpGs(-B) 

The phase shifts nf^p and df„jp are given by the integrals 



(6.28) 



2(F 



dX 

Pimp (A, Ap), 

B 
-B 

dX Pimp (A, Ap) 
dX Pimp (A, Ap) 



where pimp is the solution to 



1 - 



K 
2^ 



Pimp{X, A') — 



j^(A- A') 
2^ ' 



(6.29) 



(6.30) 



(6.31) 



Comparing Eq. (|6.26|) with Eq. (|6.6p . we conclude 

that e is the dressed energy of the single particle 



e(Ap 



(6.32) 



The parameters of the unitary transformation, which set 
the exponents of the correlation functions, are 



p 

Tr.l 



"■imp "-imp- 



(6.33) 



For a single hole, we find e = — e(A/i) > 0. The phase 
shifts for the hole case are 



■imp 



2r/'' 

%mp 



dXpimp{X,Xh), (6.34) 



-B 



dX Pimp (A, A/i) 



— oo 
oo 



+ / dXpimp (A, Xh) ■ (6.35) 

JB 

These formulas can be used to compute the exact pa- 
rameters vr^l in Eq. (I3.25P by solving the BA integral 
equation for Pimp in Eq. (|6.3ip for a given choice of parti- 
cle rapidity Ap or hole rapidity A^- In order to analyze the 



Im(Ap) 



TT 

2 







kp 



k 



FIG. 7: Momentum ranges for the three types of excited states 
considered in the BA calculations for < 1^ < 2. For |fc| < 
kp, the simplest excitation is a single hole with real rapidity 
below the Fermi boundary, \Xh\ < B. For kp < \k\ < koo, 
we add a single particle with real rapidity |Ap| > B. For 
koo < |fc| < TT, we add a particle with complex rapidity Ap = 
Re(Ap) -I- in/2 (a negative-parity one-string). 



results for correlation functions as a function of momen- 
tum k, the latter have to be chosen so that k{Xp^h) — k, 
where k{X) is the dressed momentum 

.B 

k{X)=po{X)- dA'0(A-A')pGs(A'). (6.36) 



The function k{X) is such that k{±B) = ztkp. Natu- 
rally, a single hole excitation has momentum in the range 
< \k{Xh e M)| < kp- It follows from Eq. ((OS)) that a 
single-particle excitation with real rapidity has momen- 
tum in the range kp < |fc(Ap G M)| < /cqo, where 



kF + {Tr-0[l 



2kF 
TT 



(6.37) 



A particle excitation with complex rapidity of the form 
in Eq. (|6.17p has momentum in the range kao < 
|fc(Ap = Re(Ap) + iTT/2)\ < tt. These three types of BA 
eigenstates cover the Brillouin zone, as illustrated in Fig. 
[T] The value of momentum fcoo is reached by taking 
Re(Ap) oo. From Eqs. ([QT]) and ([OT]) . we can 
see that this is the point where the energy of the single- 
particle excitation equals the absolute value of the chem- 
ical potential 



e(fcoo) = eo(A ^ oo) = |^|. 



(6.38) 



For < y < 2, we have 7t/2 < koo < tt — kp. 

Since e in Eq. (|6.32p varies continuously with Xp^h, 
the velocity u in Eq. (j3.14p is fixed by linearizing the 
dispersion of the particle excitation around Xp,h 



dk 



e'jXp.h) 
k'iXp.h)' 



(6.39) 



In Appendix ISl we show that rif^p and df^^p can also 
be expressed in terms of the shift function as done in Ref. 
[l6l . In the low-energy limit where the momentum of the 
particle or hole approaches kp < 7r/2, i.e. Xp^h B < 
oo, we use formulas (IA23P and (|A24p in the appendix and 
obtain 



p,h 



Tr 



p.h 

Tl 



±{2Vk~1 - K), 



±iK-l). 



(6.40) 
(6.41) 
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This implies that the exponent for the long-time decay 
of Gpj,{k,t) in Eqs. (jX^ and (jX^ for k ^ kp is 

Kfc^fcF<V2) = l-VX--^ + ^ + |. (6.42) 



Thus the exponent of the single-particle singularity of 
the spectral function in the low-energy limit assumes the 
value 



v{k 



< n/2) = ^+^-1^-^. (6.43) 



^/K 2K 2 



This low-energy exponent, expressed in terms of the Lut- 
tinger parameter, is actually universal (see Ref. [T3 ). On 
the other hand, we showed in Ref. 15 that in the limit 
of half-filling, kp — > 7r/2, B oo and |fc(Ap^;i)| ^ kp, we 
have 



p.h 



±{K-1) 



imp 



0, 



(6.44) 



independent of particle or hole momentum. In this limit, 
the region of validity for the particle excitations with 
real rapidities shrinks to zero as loo — > N/2 = L/A and 
koo ^ kp = 7r/2. The single hole and negative parity 
one-string are still allowed (loo ^ L/A). As a result, in 
the limit of half-filling. 



K 



1 
2^ 



(6.45) 



For comparison, the Luttinger liquid result for the par- 
ticle or hole Green's function (the Luttinger model is 
particle-hole symmetric due to the linearization of the 
dispersion) is^i^ 



G(x,t) 



-ikpx 



2'k(x — vt + irj) 2Tr{x + vt — irj) 

^2 



x'^ — {vt — irjy 



(6.46) 



Taking the Fourier transform to momentum space, we 
get the long-time decay 



G{k w kp,t) - r 



(6.47) 



Note that the exact exponent for the decay of 
Gp.h {k ~ kp,t) in Eq. (|6.42p is smaller than the Lut- 
tinger liquid exponent vu. The Luttinger liquid result 
for the singularity of the spectral function on the single- 
particle energy 



A{k,ujK v6k) - Sk''"/^{uj - vdky^+''"/'^ 



(6.48) 



where (5fc = fc—fci?. Interestingly, according to Eq. (|6.45|) . 
in the limit of half-filling the exact exponent v of the 
long-time decay of Gp(fc, i) (for |fc| > kp) or Gh{k,t) (for 
I A; I < kp) becomes independent of k (for any k in the 
entire Brillouin zone) and coincides with the Luttinger 
liquid exponent vu . However, the exact exponent for 



the single-particle singularity of A{k,uj) still differs from 
the Luttinger liquid result (note the factor of 1/2 in Eq. 
([OSll compared to Eqs. ([^3T|l . (|3:32l) and ([335]) ). The 
difference stems from the fact that in the Luttinger model 
the singularity of the spectral function at cj vSk is 
controlled by the singularity of the Green's function in 
the vicinity of the light cone x « — wt, rather than by the 
exponent for t ^ \x\. 

That the half-filling exponent can be expressed solely 
in terms of the Luttinger parameter is a peculiarity of 
the integrable model. The exception is the model with 
V — 2, which is equivalent to the Heisenberg spin chain. 
If additional interactions which break integrability still 
respect SU(2) symmetry, the exponent = i^u = l/A for 
K ^ 1/2 is universal.— 



B. One- two-string one- hole excitations 

The Bethe equations in Eq. (|6.7p admit solutions with 
a pair of complex rapidities of the form 



As=w + </2, X*^w~iC/2, wG 



(6.49) 



An eigenstate with two particles described by such pair 
of rapidities is called a two-stringi^ The fact that 
gje(As-Aj _ g guarantees that the wave function of the 
two-string is normalizable and describes a two-particle 
bound state. 

From Eqs. (|6.8p and (|6.1ip . we have that the bare 
momentum and energy of the two-string are 



Voiw) =Po(As) +Po(A*) = ilog 



sinh(i^ + w) 
sinh(iC — w) 



£o{w) ^ eoiXs) + eo{X:) 



2 cos ( sin C 
cosh^ w — cos^ C 



(6.50) 



2^1. 



,;.5i) 

Note that the bare momentum of the two-string is re- 
stricted to the interval 



iVoiw) - 7r| mod 27r < tt - 2C. 



(6.52) 



If the two-string is added to the ground state of N 
fermions, the renormalized momentum and energy must 
include the backflow of the ground state. In analogy 
with Eqs. (|6.27p and (|6.36p . the dressed momentum and 
energy of the two-string excitation are 

r + B 

V{w) = Voiw) - / dA'e(u.- A')p(A'), (6.53) 



£{w) = £o{w)+ / /C(w; - A')e(A'), (6.54) 

J-B 27r 



where 



e(w) = eiw + iC/2) + 9{w-iC/2), (6.55) 
JCiw) = K{w + iC/2) + K{w~iC/2). (6.56) 
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One can then verify that the dressed momentum is re- 
stricted to 



iViw) - 7r| mod 27r < Abs, 



where 



A, 



7r-P(-oo) = (7r-2C) 



(6.57) 



(6.58) 



This fixes the momentum range for the bound state band 
alluded to in Sec. IIVI The minimum value of momentum 
k for which a one-two-string one-hole excitation exists is 



kh — ki 



2kF 
2C 1 - — 

TT 



(6.59) 



Moreover, it is easy to show that the energy of the bound 
state at the edge of the bound state band is 



£{V = n± Abs) = £{w ±oo) = 2\n\ 



(6.60) 



In the limit of half-filling, kp 7r/2, the bound state 
band shrinks to zero {Abs — > 0) for arbitrary values of 
the interaction strength and the region of validity of the 
power-law singularity vanishes. 

The phase shifts at the Fermi boundaries due to the 
creation of a two-string can be calculated by the methods 
described in the previous section. Defining ft„ip(A, u) by 
the integral equation 



IC{X-w) 
2^ ' 



(6.61) 



we can write the phase shifts that enter the finite size 
spectrum as 



2d'. 



dXgimp{X,wo), (6.62) 





r+B 


bs _ 




'imp 


/ ' 


lbs 




'imp 


J -—oo 

Jb 



d\Qirnp{X,Wo), (6.63) 



where wq is chosen such that V{wf)) = Pq for a bound 
state with center-of-niass momentum Pq. By correspon- 
dence with the finite size spectrum of the field theory, the 
parameters of the unitary transformation that determine 
the exponents in Eq. (|4.36p are given by 



r^bs 



j6s 



imp ^imp' 



(6.64) 



The velocity of the bound state is obtained by lineariz- 
ing the dispersion about wq 



Ubs = 



(6.65) 



VII. DENSITY MATRIX RENORMALIZATION 
GROUP RESULTS FOR FERMION GREEN'S 
FUNCTION 

A. Method 

We used real time DMRG methods (tDMRG) to cal- 
culate the spectral functions for this system. The typical 
size of the systems studied was L = 400. The main de- 
tails of the techniques used were described briefiy in Ref. 
il5i, and in considerable detail in Ref. i4Q . Here we sum- 
marize the method, and also discuss a finite size issue 
which did not come up in the earlier work, and how we 
treat it. 

After finding the ground state with ground state 
DMRG, an operator ■0 or -ip^ was applied to a site in the 
center, and the resulting solution to the time-dependent 
Schrodinger equation was used to calculate a space-time 
dependent hole Green's function 



G;,(a;,i) = (0]+,(i)0^.(O)) 
or particle Green's function 

Gp(a;,i) = (^,+.(t)V'l(0)) 



(7.1) 



(7.2) 



for t > out to a particular time T ^ 20 — 30. The 
hole and particle Green's functions each require a sepa- 
rate simulation. The result was Fourier transformed in 
space. For each desired wave vector fc, the resulting time- 
dependent functions Gp^h{k,t) as defined in Eqs. p.22p 
and (j3.33p were either fit to an asymptotic long-time form 
and extrapolated to very long times using it, or extrap- 
olated to long times using a very general method called 
linear prediction. In both fast Fourier trans- 

form was used to obtain high resolution spectra. The 
asymptotic form also gave directly exponents character- 
izing the frequency singularities. Typically, we speci- 
fied a truncation error of 10~* to set the number of 
states kept per block m, with an additional constraint 
m < 1500. These parameters were varied to determine 
errors in Gpji{x, t), and the typical errors for the parame- 
ters used were 10""* — 10~^. The functions Gp^h{x, t = Q) 
decay as a power law in x, potentially giving large finite 
size effects. However, these tails in x only appear in the 
real part of Gp^h, and it is easy to reconstruct the entire 
spectrum from only the imaginary part when the particle 
and hole Green's functions are treated separately. The 
support of the imaginary part spreads out from the cen- 
ter site with a speed given by the maximum velocity Vm 
as a function of k of the hole or particle; finite size effects 
are small as long as Vm,T ^ L/2. The long tails in x 
are associated with w « 0; the determination of spectra 
using only the imaginary part generates an odd function 
in frequency, eliminating the a; = contribution. 

Another finite size effect comes from the Friedel os- 
cillations induced by the boundaries of an open system 
in the ground state. In Fig. [SJ the solid dots represent 
the ground state density on each site {n{x)) of an open 
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FIG. 8: (Color online) Local density {n{x)) in the ground 
state for an average density of n = 0.45 and interaction 
strength V = 1, near the the center of a system with sys- 
tem size L = 400. Black dots: open boundary conditions; 
solid red line: smooth boundary conditions. 




0.5 - 



400 



FIG. 9: Smoothing function s{x) used in the DMRG runs 
with smooth boundary conditions. 



L = 400 system with a desired average density of 0.45. 
The open boundaries induce oscillations in the density 
which decay slowly away from the ends. The density is 
directly related to Ghix = 0,t = 0), and these oscilla- 
tions introduce a finite size error in this system of order 
1%. There are several ways one could substantially re- 
duce the finite size effects. The simplest (conceptually) 
would be to use periodic boundary conditions, which can 
now be treated within DMRG with only a slight compu- 
tational penalty.'*^ Because our most efficient computer 



program for real time evolution does not incorporate the 
new periodic algorithm, we have chosen to use smooth 
boundary conditions (BCs)4^ For smooth BCs, one in- 
troduces a chemical potential to set the density in the 
center of the system, but one still works with a fixed 
number of particles. The parameters of the Hamiltonian 
and the chemical potential are held fixed at bulk values 
in a central region of the system, but near the edges they 
are scaled down to zero smoothly; see Fig. |9l The basic 
smoothing function by which the parameters were mul- 
tiplied (after scaling and shifting) was 



s(x) 



^|l + sin |cos(7ra;) |, (7.3) 



which falls smoothly from 1 to over < a; < 1, and 
has three vanishing derivatives at the endpoints. The 
smoothing regions were one fourth the lattice length at 
each end. The smooth boundaries allow the central bulk 
region to minimize its energy, largely eliminating the 
Friedel oscillations, at the expense of large oscillations 
near the low-energy edges. The chemical potential allows 
one to fix the bulk density to any value, even irrational. 
One sets the overall number of particles to the integer 
giving the closest overall density, and the system adjusts 
the density at the edges to give a density in the center 
determined by the chemical potential. (This means that 
a number of ground-state-only runs must be done to de- 
termine the correct chemical potential for the given den- 
sity; these contribute little to the overall computational 
cost.) The smooth connection between the boundary re- 
gion and the bulk is essential to avoid artifacts there. 
Fig. [5] shows the resulting ground state density, with ac- 
ceptably small oscillations in the central region. The real 
time evolution is performed using the smoothed Hamilto- 
nian. The evolution should be stopped before the signal 
initiated by ip or tp^ in the center reaches the edge of the 
bulk region. We have checked how well the smooth BCs 
perform in determining the correct thermodynamic limit 
Green's function in the case V = 0, where one can 
obtain numerically exact results from diagonalizing and 
manipulating L x L matrices. We find finite size errors 
of 10^'^ — 10~^ for the size systems we use. 

The smooth BCs duplicate the behavior of a very large 
system in another respect: the entanglement entropy, 
which governs the number of states kept per block m for 
a given error, is significantly larger than with open BCs 
for a given L. We have not studied this interesting effect 
in detail. We may obtain a simple understanding of this 
effect by considering two weakly coupled spins at the two 
ends of a Heisenberg antiferromagnetic spin chain of even 
length. If these end spins are sufficiently weakly coupled 
then they will form a singlet together in the ground state, 
whereas all other spins go into a complicated "resonating 
valence bond" state. This long range singlet increases the 
entanglement of the left and right halves of the system. 
Related entanglement phenomena in the case of a single 
weakly coupled spin were studied in Ref. IH. We were 
able to increase the number of states kept sufficiently to 
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TABLE I: Parameters for the model with V = 1 and n — 
0.2, obtained either from formulas given in Sec. IVII or from 
the numerical solution of the BA equations for finite Fermi 
boundary B. 



K 



Po - kf 



-2.542 1.343 0.873 0.628 1.885 1.885 2.714 



obtain accurate spectral functions despite the larger en- 
tanglement. An intuitive picture explaining this effect 
is that the smooth wave function might be written as 
a superposition of a non-translationally invariant wave 
function, such as shown by the circles in Fig. [51 with 
translations of itself, sufficient to eliminate the static os- 
cillations. If this idea is roughly right, then one might 
obtain accurate spectra by using one open BC lattice, 
but averaging over the starting point where ip or -0^ is 
initially applied, say with a Gaussian envelope of width 
5 — 10, to smooth out the effect of the oscillations. We 
have tried this approach as another check on the smooth 
BC results. It requires at least 5-10 runs, one for each 
starting point, but each has low entanglement. The re- 
sults were quite satisfactory, giving good agreement with 
the smooth BC results. 



B. Results 

We used the BA integral equations in Sec. |Vl]to nu- 
merically evaluate the dressed energies e(A) and S{w) of 
the elementary excitations for the model with V = 1 
and n = 0.2. The values of some important parameters 
are presented in Table HI Note that for y = 1 we have 
kb = fcoo, thus the range of momentum where the bound 
state contribution exists coincides with the range where 
the single-particle excitation is described by a negative- 
parity one-string. The result for the dispersion of the 
single-particle excitations and the lower threshold of the 
one-two-string one-hole continuum is shown in Fig. 1101 
We find that, for n = 0.2 and V = 1, the velocity of 
the two-string equals the renormalized Fermi velocity at 
Pq « 3.342 < TT + Abs. As a result, the nature of the 
bound state singularity changes at fc = Pq* — kp ^ 2.714 
(see Fig. [TU]). 

Using the tDMRG method, we calculated both par- 
ticle and hole Green's function for the values n — 0.2 
and y = 1 and for times out to T = 37. The tDMRG 
results show that the particle Green's function Gp{k,t) 
becomes 0(1) for |fc| > kp- Due to interactions, Gp{k,t) 
is nonzero for momenta below the Fermi surface, but its 
amplitude is much smaller. Similarly, the hole Green's 
function is 0(1) for |fc| < kp- We restrict our analysis to 
the ranges \k\ > kp for the particle Green's function and 
|fc| < kp for the hole Green's function. As an example, 
Fig. [TT] shows the tDMRG result for the particle Green's 
function Gp{k,t) for k — 0.357r. This is in the fast par- 
ticle regime u > v. The numerical result is compared 
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FIG. 10: (Color online) Exact single particle dispersion and 
lower threshold of the one-two-string one-hole continuum for 
interaction strength V — 1 and fermion density n = 0.2. The 
curve uj = e{k) is divided into the hole region (fc < kp) and 
two particle regions (particle excitations with real rapidities 
for kp < k < koo and complex rapidities for koo < k < tt). 
The dashed line is a cosine function fitted to the hole region 
of the spectrum, showing that the renormalized dispersion 
deviates from the cos fc dependence. The one-two-string one- 
hole continuum is defined for kb < k < n (see text). For 
kb < k < Pq — kp , the lower threshold is defined by a hole at 
kp and a two-string with momentum k + kp. For Pq — kp < 
k < TV, the lower threshold is defined by a deep hole and a 
two-string with the same velocity. 



with the field theory result of Eq. (|3.28p . which pre- 
dicts an asymptotic behavior with oscillations at a single 
frequency and pure power-law decay. The BA fixes the 
frequency, phase and exponent of the decay of Gp{k,t). 
Such very good agreement of Eq. (|3.28p with the tDMRG 
result is typical for k values in the regime u > v. This 
lends support to the conjecture that the decay rate 1/r 
discussed in Sec. IV Al is exactly zero for the integrable 
model. 

Similar results with a pure power-law decay are ob- 
tained for the hole Green's function Gh{k, t). In this case 
we expect 1/t to vanish due to kinematic constraints, 
since the single hole is at the threshold of the multipar- 
ticle continuum. (The exact edges of the support can be 
determined from the dispersion in Fig. 1 101 and are similar 
to the curves in Fig. [2^.) 

For contrast, we show in Fig. [T2|the tDMRG results 
for the particle Green's function for k = O.SStt. The data 
is typical of that for values of k near tt and shows the 
presence of two dominant frequencies in the long-time 
behavior. 

By utilizing the linear extrapolation method described 
in Ref. 40, we computed the Fourier transform of the real 
time tDMRG data to produce line shapes for the spec- 
tral function without any analytical input. We confirm 
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FIG. 11: (Color online) Real part of the particle Green's func- 
tion Gp{k — 0.357r,t) for V = 1 and n = 0.2. Red dots are 
numerical tDMRG results and solid line (FT+BA) is the field 
theory result of Eq. (|3.28p with parameters computed from 
the Bethe ansatz: frequency e^^ — 0.8105, exponent = 
0.0097 and phase ly^^ - = 0.0076. Only the amplitude 
has been fitted. If we fit the DMRG results for i > 5 with 
four free parameters, we find the best fit for g^*^^"^ = 0.8104, 

^DMRG ^ o^og ^DMRG _ ^DMRG ^ qqqj^ 




FIG. 12: (Color online) tDMRG result for the real part of the 
particle Green's function for V — 1, n — 0.2 and momentum 
k = O.SStt. The line is a guide to the eye. 



the simple behavior of the hole spectral function, with 
a single sharp peak on the single-hole energy. We thus 
focus on the particle part of the spectrum, which shows 
a more interesting behavior. The results for k > 0.357r 
are shown in Fig. [131 The energies of the peaks in the 
tDMRG results agree with the energies predicted by the 
BA solution shown in Fig. [TOl The single-particle peaks 
remain sharp for fairly large values of k — kp which are 
in the regime u > v {k < 2.31 according to the BA dis- 
persion in Fig. [TU)l . On the other hand, the broadening 
of the single-particle peak is apparent for the larger val- 
ues of fc, near fc = tt. This happens in the slow particle 
regime u < v, where we expect a nonzero decay rate due 
to the process discussed in Sec. IV Bl 
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FIG. 13: (Color online) Particle spectral function Ap{k,uj) 
for V = 1, n = 0.2 and k/iv = 0.35, 0.4, 0.45, 1, obtained 
by Fourier transforming the tDMRG data using the linear 
extrapolation method. The curves for different values of k 
are shifted vertically {k increases from bottom to top). The 
horizontal dotted line (fc ~ 2.31) separates the regime u > v 
from the regime u < v. For kp < k < kb, there is only 
one peak at the energy of the single-particle excitation. For 
kb < k < 77, there appears a second peak, which we interpret 
as due to a two-particle bound state and a free hole. The 
energy of the second peak is non-monotonic in k and has 
a minimum at k = n — kp- The dashed lines are the BA 
predictions for the dispersion relations in Fig. 1101 



The second peak associated with the one-bound-state 
one-hole excitation shows up for momentum k > kb, as 
predicted by the BA. For k near kb, the second peak is 
small and broad and hardly visible in Fig. [T31 This may 
be explained qualitatively by arguing that, when the ve- 
locity of the two-particle bound state is smaller than the 
Fermi velocity, the decay of the bound state by scatter- 
ing of one low-energy particle-hole pair is kinematically 
possible (for reasons analogous to the arguments given 
in Sec. IV Bp . so the edge singularity may actually be 
broadened. 

As k increases towards k — spectral weight is trans- 
ferred from the single-particle peak to the bound-state 
hole continuum. The second peak also becomes sharper 
with increasing k, in accord with a stronger power-law 
singularity in the regime where the lower threshold is de- 
fined by a bound state and a deep hole with the same 
velocity (see Sec. lIVp . 

We obtain high resolution spectral functions by using 
the asymptotic behavior predicted by the field theory to 
extend the tDMRG data to arbitrarily large t before com- 
puting the Fourier transform. The result of this combi- 
nation of methods is illustrated in Fig. [14] for the particle 
spectral function in the fast particle regime, u > v. In 
Fig. [iniwe show the hole spectral function for a different 
value of density, n = 0.45. (For n = 0.2, the maxi- 
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FIG. 14: (Color online) Spectral function in the vicinity of 
the single-particle energy for 1/ = 1, n = 0.2 and k = O.SStt, 
obtained by extending the tDMRG data shown in Fig. Illl to 
long times using the field theory formulas and then taking the 
Fourier transform. 




FIG. 15: (Color online) Spectral function in the vicinity of 
the single-hole energy for V — 1, n — 0.45 and k = O.Olvr, 
obtained by extending the tDMRG data as in Fig. 1141 



mum energy of a hole is small, so the tDMRG data with 
t < 40 is not representative of the long time behavior.) 
While the spectral function exhibits a two-sided power- 
law singularity on the single-particle energy, it vanishes 
below the single-hole energy. These line shapes represent 
a high-energy extrapolation of the low-energy result of 
Khodas et al.,1^ except for the absence of broadening of 
the single-particle peak. We are not able to access di- 
rectly the regime k kp because in the low energy limit 
the region of validity of the power law in Eq. p.3ip be- 
comes extremely small. This means that it would be nec- 
essary to go to extremely large values of t in the tDMRG 
calculations in order to observe the asymptotic behavior 
of Gp{k,t). The contribution of the convergent edge sin- 
gularities to the long-time decay of Gp{k, t) is also rather 
small, so we have not attempted to extract the behavior 
near the multiparticle thresholds where A{k, uj) vanishes. 
We have also calculated the particle Green's function 



for the same V = 1 and n = 0.45, which is closer to 
half-filling. In this case both hole and particle spectral 
function show a single sharp peak a,t u! — e{k). We have 
not observed a second peak in the particle spectral func- 
tion for k near tt. This is in agreement with the prediction 
that the contribution of the one-bound-state one-hole ex- 
citation disappears as we approach half-filling (see Sec. 

MM- 

Assuming a single frequency is able to describe the 
long-time decay of particle and hole Green's function, 
we fitted the data to Eq. (|5.2ip in order to extract a 
possible nonzero decay rate. Eq. (|5.2ip is valid for k 
values in the regime u < v, in which a single parameter 
fixes both phase and exponent. The results are shown in 
Fig. \M We first note that fitting the tDMRG data up to 
a maximum time T = 20 is not reliable near two points. 
The first point is k — kp, where the frequency of the time 
oscillation vanishes. The second point is fc ~ 1.897, where 
u~v according to the BA. Near this point, the time scale 
for reaching the asymptotic decay also diverges (see Eq. 
()3.27|) ). These two points are indicated by vertical lines 
in Fig. [TH The agreement between the tDMRG and BA 
exponents is expected to be best near k = for the hole 
case and near k = ir for the particle case. 

For the hole Green's function, we find very good agree- 
ment between the numerical exponent and the one pre- 
dicted by the BA. Moreover, in this case the numerical 
decay rate is negligible (< 10^^ for all values of A; < kp). 
For the particle Green's function, the agreement between 
BA and tDMRG exponents is not as good as for the hole 
case. This is expected given that we are fitting with two 
small parameters, the power law decay exponent and the 
decay rate. We have checked that fitting the tDMRG 
results with the exponents constrained by the BA works 
as well as fitting with unconstrained exponents. Thus 
the oscillations in the tDMRG exponents in Fig. [16] are 
believed to be due to the error in the fitting. Nonethe- 
less, the tDMRG results clearly show that the decay rate 
1/t is nonzero for a high-energy particle in the regime 
u < V. The results also suggest that 1/t decreases as k 
approaches the value where u = v, possibly vanishing at 
the lower threshold. We note that the behavior of the BA 
exponent is qualitatively similar to the expected from the 
weak coupling expressions in Eq. (I3.17p . 

Finally, we point out that the long-time behavior of 
the fermion Green's function in real space is dominated 
by a saddle point contribution^^ The latter corresponds 
to a hole at the bottom of the band, fc = 0, in the case 
of the hole Green's function, or to a particle at the top 
of the band, fc = tt, in the case of the particle Green's 
function. The high-energy excitations near these points 
have parabolic dispersion e(fc) « e + k'^/2m, where m < 
is the renormalized mass. Thus the propagator for the 
decoupled d particle in Eq. (|3.26p reads 



(:r,t) 



2nit 



(7.4) 



Therefore, in the noninteracting case, this contribution 
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FIG. 16: (Color online) Exponent of particle (for k > kp) 
or hole (for k < kp) Green's function for n = 0.45 and V = 
1. The solid lines represent the exponents extracted from 
tDMRG by fitting the data in the time interval 10 < t < 20 
to Eq. (|5.2ip . The dotted line is the decay rate of the single 
particle extracted from the fitting. The dashed lines represent 
the exponent calculated numerically from the BA. The dot at 
k = kp indicates the universal low-energy exponent in Eq. 
(I6.42p . The vertical lines indicate the values of momentum 
k ^ kp = 0.457r and k ^ 1.897, near which the tDMRG 
exponents are least reliable. 



to the particle Green's function oscillates at a high fre- 
quency 0(1) for t 3> \m\x^ and decays as l/\/t. In 
the interacting case, the frequency of the oscillation be- 
comes the renormalized energy of the single particle with 
momentum k — tt. In addition, the exponent is modified 
by the coupling to the low-energy modes. In the case 
of the particle Green's function, there is also an expo- 
nential decay associated with the decay rate of the 
particle at A: = tt when the model is below half-filling. 
The long-time behavior of the particle Green's function 
is then given by 



Gp{x, t ^ x/v, \m\x'^) 



B 



(7.5) 



Here A and B are complex amplitudes. The second term 
in Eq. ()7.5|) is the standard Luttinger liquid result, which 
does not oscillate in time. The renormalized exponent of 
the high-energy term is 



ry^ = l/2 + iy{k = 7T), 



(7.6) 



with u{k) given by Eqs. ([X^ and f^^ . This is 
smaller than the exponent of the Luttinger liquid term; 
however, for t > the high-energy term in Gp{x, t) 
decays exponentially. We note that both a particle at 
k = IT and a hole at A; = couple symmetrically to 



TABLE II; Exponents and frequencies for the fermion 
Green's function G{x = 0, t) at half-filling and for V = 
0.25,0.5,0.75,1,1.5. The parameters edmrg, Vdmrg and 
ctDMRG were obtained by fitting the DMRG results to Eq. 
(|7.1ip . The other parameters are analytical BA predictions 
and should be compared with the adjacent DMRG results. 



V 


K 


£dmrg 


eo 


Vdmrg 


'70 


0!DMRG 


l + t-H 





1 




2 




0.5 




1 


0.25 


0.926 


2.156 


2.156 


0.503 


0.503 


0.985 


1.003 


0.5 


0.861 


2.308 


2.308 


0.510 


0.511 


0.986 


1.011 


0.75 


0.803 


2.454 


2.454 


0.523 


0.524 


1.002 


1.024 


1 


0.75 


2.598 


2.598 


0.539 


0.542 


1.013 


1.042 


1.5 


0.649 


2.876 


2.876 


0.583 


0.595 


1.048 


1.095 



the low-energy modes at the two Fermi points, hence 
— 0,Tr) = VL{k — 0, 7r). Below half-filling the hole 
Green's function decays as 



Gh{x, t ^ x / V , \m\x'^) 



A' 



tVo 



B' 



(7.7) 



with 770 — 1/2-1- v^k — 0) and no exponential decay of 
the high-energy contribution. 

At half-filling, particle and hole Green's functions are 
equivalent. In this case, we have analytical expressions 
for the frequencies and exponents^- 



^ arccos(l//2) 
r/^ = 1/2 + i^u, 



(7.8) 
(7.9) 



where i^n = [K + K ^ — 2)/2 with the Luttinger 
parameter— 



K = 



2[7r-arccos(y/2)]' 



(7.10) 



Moreover, the decay rate 1/t^ vanishes at half- filling. 

We have calculated the fermion Green's function 
G{x = 0, t) at half-filling and for various values of in- 
teraction strength V using tDMRG. We used simulated 
annealing to fit the numerical results to the formula 



Gix,t) 



pi DMRG 



B 



fO^DMRG 



(7.11) 



The results are shown in Table ITTl We find good agree- 
ment with the analytical BA predictions of Eqs. (|7.8p 
and (|7?gi) for all values of V. 



VIII. DYNAMICAL STRUCTURE FACTOR 

The methods used in this work can be generally ap- 
plied to the study of singularities of dynamical correla- 
tion functions. The basic idea is to identify the excita- 
tions which define the thresholds of the spectrum in the 



24 



noninteracting limit and then ask how interactions be- 
tween the particles in the final state affect the behavior 
of the dynamical function near the threshold. This can 
be answered with the help of effective models such as 
Eqs. (|3.14p . One quantity of interest is the dynamical 
structure factor 

S{q,u;) = ^e-'^O-^- ) / dt e^-\n,{t)n,, (0)) . 

3,3' 

The dynamical structure factor for the Galilean invari- 
ant model was studied in Ref. [l^- In the lattice model, 
S{q,ijj) is equivalent to the longitudinal dynamical spin 
structure factor of spin- 1/2 chains.— A particular inter- 
esting case from the point of view of experiments is the 
model with V = 2^ which is equivalent to the Heisenberg 
spin chain 

L 

H = jY,SrS,+i. (8.2) 

3 = 1 

In this section we discuss qualitatively the implications 
of our results on the spectral function for the line shape 
of S{q, uj) away from half-filling. 

The excitations which give the dominant contribu- 
tion to the dynamical structure factor were identified by 
Miiller et al4^i^ In the noninteracting case, the region 
where S{q, w) is nonzero is given by the particle-hole con- 
tinuum shown in Fig. [17] for n = 0.4. For small momen- 
tum, q <C TT — 2kp, S{q,uj) is similar to the result for 
the Galilean invariant model »2£ Instead of discussing all 
possible cases, here we focus on the momentum range 
IT — 2kp < q < 2kp. In this case there are three impor- 
tant thresholds. Let us denote by k the momentum of 
the hole, so the momentum of the particle is k + q. The 
lower threshold uobiq) of the particle-hole continuum is 
defined by a deep hole excitation composed of a particle 
with momentum kp and a hole with k — kp — q. The 
upper threshold i^uil) is defined by the particle- hole ex- 
citation which is symmetric about 7r/2, i.e., k — {■n—q)/2. 
There is yet a third threshold uj]\i{q), situated between 
LOL{q) and LOu{q), defined by a high-energy particle with 
momentum k + q ^ kp + q and a hole at kp. 

The exact S{q, lu) for the noninteracting model in the 
range tt — 2kp < q < 2kp ia^ 

^ 9{^u{q) - - iodq)) + - ^Af (g))] 

yjuj'fjiq) - Cj2 

(8.3) 

According to Eq. (|8.3p . S{q,Lo) has step discontinuities 
both at uJhiq) and ujM{q)- The spectral weight jumps by 
a factor of two at uj]\i{q) because for lo > LO]\i{q) there 
are two choices of hole momentum k corresponding to the 
same energy. For ^^(q) < < cjmC?), there is only one 
choice of k for each uj. At the upper threshold LUij{q), 
S{q,uj) has a square root divergence which stems from 
the divergent joint density of states of a particle and a 
hole with the same velocity. 




FIG. 17: Particle-hole continuum for the noninteracting 
model with n — 0.4. The dynamical structure factor S{q,u!) 
is nonzero in the shaded area. 



In the limit of half- filling, kp — > 7r/2, the deep hole 
excitation and the high-energy particle excitation for the 
same q become degenerate. The intermediate threshold 
merges with the lower threshold and one recovers a single 
step function at WL(q)4^ 

By analogy with the ansatz for the Heisenberg chain 
at zero magnetic field, Miiller et al.^** proposed that in 
the presence of repulsive interactions the two step dis- 
continuities should become divergent edge singularities. 
This was argued to explain the double peak structure 
observed in the dynamical structure factor of spin chain 
compounds at finite field;^^ 

The picture of a line shape with two peaks at wl(<z) 
and LOM{q) is supported by the field theory approach 
to the edge singularities of S{q,uj). The divergence at 
the lower threshold is easily understood as an x-ray edge 
singularity for the deep hole excitationiiS For the upper 
threshold, it is known that, at half-filling, repulsive inter- 
actions turn the square-root divergence into a universal 
square-root cusp.^^ The reason is the resonant scattering 
between the high-energy particle and hole with the same 
velocity, which is analogous to the problem of Wannier 
excitons in semiconductors.— Away from half-filling, the 
square-root divergence still becomes a convergent singu- 
larity. However, the behavior near u>u{q) is not exactly a 
square-root cusp because in the absence of particle-hole 
symmetry the high-energy particle-hole pair does not de- 
couple from the Fermi surface modes»iS, 

The behavior near LJM{q) is controlled by the x-ray 
edge singularity of a high energy particle excitation with 
u < V. Using the methods of Ref. [13, it is easy to show 
that S{q, oj) acquires a one-sided divergent power-law sin- 
gularity above uJMiq)- The exponent of this singularity 
is different from the one at the lower threshold. Using 
the weak coupling expression for the phase shifts, we can 
show that the exponent at LOM{q) to first order in V is 




given by 

/ N ^(1 ~ COS 9) 

7r[sin kp — sin(fci? + q)\ 

On the other hand, the exponent at ojLiq) to 0{V) is 

. ^ y(l -cosg) 

7r[sinA;i? — sm(kp — q)\ 

For g > TT — 2fci?, we have 

t^Aiiq) > Ml('7)- (8-6) 

However, based on the results of Sec. I VIII we expect 
that below half-filling the high-energy particle has a finite 
decay rate 1/r when u < v. Therefore, we predict that 
S{q,uj) has a rounded peak rather than a divergence at 
(g). A schematic picture of this line shape is shown 
in Fig. \M 

In the limit kp — > 7r/2, the peak at ujM{q) approaches 
the divergent singularity at ujL{q). At the same time, 
liM becomes equal to ^l, as required by particle-hole 
symmetry. In addition, at half-filling the decay rate of 
the high-energy particle vanishes (see Sec. IV Bp . As 
pointed out in Ref. Tg*, the range of validity of the ex- 
ponent /ii shrinks to zero as kp — > 7r/2. This was ar- 
gued to lead to a discontinuity in the exponent, such 
that fip^kp — 7r/2) ^ ^^L{kp — > 7r/2). However, what 
happens is not a crossover between two power-law sin- 
gularities, as in the case of the spectral function,— but 
rather a collapse between two divergent singularities with 
the same exponent. Our scenario predicts that S{q,uj) 
varies smoothly as kp ^ it/2. At kp = Tr/2, one re- 
covers the line shape with a momentum independent ex- 
ponent at the lower threshold, = 1 — K.i^ Indeed, 
the results of Ref. [is' show that the alternative expo- 
nent fiLikp = 7r/2) proposed in Ref. [l^ is inconsistent 
with universal relations for the phase shifts in the low- 
energy limit and with the SU(2) symmetry of the model 
for V = 2. The correct exponent at half- filling is the one 
derived in Ref. [Tsl . 

The line shape of S{q,uj) shown in Fig. [TH] can be 
tested by tDMRG calculations for the spin-spin correla- 
tion function of the XXZ model at finite magnetic field 
or by direct computation of form factors from the BA, 
including negative-parity one-string excitations i'^^'^^ 

Here we have only discussed the nature of the thresh- 
olds defined by free particles and holes. We note that 
a recent BA work has shown that two-string excitations 
(two-particle bound states) give an important contribu- 
tion to the transverse dynamical structure factor of the 
spin-1/2 chain at finite magnetic fieldi^ 

IX. CONCLUSION 

We have studied the spectral function of spinless 
fermions on a ID lattice using a combination of field the- 
ory methods, Bethe ansatz and time-dependent DMRG. 



FIG. 18: Schematic illustration of the line shape of S{q,uj) 
for TT — 2fcF < q < 2kF and small V. The most prominent 
features are the divergent singularity at the lower threshold 
and the rounded peak at the intermediate threshold. 



We showed that x-ray edge type effective models are 
able to explain the long-time decay of the single-particle 
Green's function and the singularities of the spectral 
function calculated by tDMRG. 

The results can be summarized as follows. In the lat- 
tice model, the detailed line shape of the spectral function 
depends strongly on momentum, density and interaction 
strength. Kinematically, the support of the spectral func- 
tion has finite-energy lower thresholds for general mo- 
mentum k only if kp /tt is a rational number, likp/ir is ir- 
rational, so that arbitrary Umklapp scattering processes 
are incommensurate, the spectral function is nonzero at 
all energies. In any case, a general approximate picture 
for the line shape can be obtained by focusing on the 
possible divergent singularities on the single-particle en- 
ergy or near multiparticle thresholds which involve bound 
states. 

Away from half-filling, the hole and particle contribu- 
tions to the spectral function are remarkably different. 
For kp < 7r/2, the hole spectral function exhibits a sin- 
gle divergent singularity above the single-hole energy. By 
contrast, the particle spectral function can exhibit one or 
two peaks. The single-particle peak is always present. In 
the regime where the velocity of the high-energy particle 
is larger than the Fermi velocity (u > v), the field the- 
ory predicts a divergence from both sides of the single- 
particle energy. For a generic model, this singularity is 
expected to be replaced by a broadened peak due to de- 
cay processes involving three-body collisions. For the 
integrable model we have found no compelling numerical 
evidence of broadening of the single-particle peak in this 
regime, even at fairly high energies, suggesting an exact 
power-law singularity. 

However, at higher values of momentum, for which the 
velocity of the high energy particle is smaller than the 
Fermi velocity, we found that the single-particle peak 
has a nonzero broadening. We interpret this in terms 
of a decay rate l/r which is second order in the interac- 
tion strength. The broadening is expected to occur for 
integrable as well as non- integrable models for u < v. 
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In real time, the corresponding particle Green's function 
decays exponentially. 

The second peak in the particle spectral function ap- 
pears at high energies, above the single-particle energy, 
and is more pronounced for low densities and momen- 
tum near tt. This peak has a non-monotonic dispersion 
relation, with an energy minimum at k — t: — kp. The 
nature of the second peak can be understood as the sin- 
gularity at the lower threshold of the continuum defined 
by a two-particle antibound state and a free hole. 

The results of this work could be tested by measuring 
the spectral function of fully spin-polarized fermions in 
optical lattices. The most robust effect should be the 
observation of the bound state peak. In a photoemission 
experiment, where one probes the hole contribution to 
the spectral function, the bound state peak should be 
visible at high densities (above half-filling) and for strong 
interactions. In this case, it can be interpreted as the 
process in which the outcoupled atom leaves behind a 
hole plus a particle-hole excitation. While the particle 
propagates freely, the two holes form a stable repulsively 
bound pair in the p-channel. 
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APPENDIX A: FINITE SIZE SPECTRUM 

In order to derive the finite size spectrum from the BA 
equations, we start by expanding the sum in Eq. (|6.19p 
for large L using the Euler-Maclaurin formula^ 



P(A) 



27r 



24L2 



K'{\ - B 

P{B+) 



K{X~X')p{X' 
K'{\ 



^P.h (A) 



27ri 



P{B-] 



(Al) 



where we introduced the Fermi boundaries B± — A(/±). 
We organize the solution to Eq. (jAip by orders of 1/L. 
We expand p(A) up to 0{1/L?') in the form 



p(A) = poo(A|B+,B_)± 



Pimp (A, \p^h\B+, 



B) 



PiiX\B+,B_) p,{-X\-B_,-B^ 
24:L^p{B+) 24L2p(S-) 



where the plus sign is for pimpi^, ^p\B+., -B_), in the case 
of a particle, and the minus sign for Pimpi^, A/i|i3-|-, B-), 
in the case of a hole. Substituting Eq. (|A2|) into Eq. 
(|A1|) . we find the integral equations for the terms in the 
expansion of p{X) 



J 



Poo(A|S+,S_) 
PimpiX, X'\B+,B^) 
p,{X\B+,B^) 



P'oW 



W+ r^^if(A-A')Poo(A'|i?+,B_), 



K{X-X) 
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^+ dX" 

— K{X-X")p.„np{X",X'\B+,B^), 
I Zn 



Ztt Jg_ In 



(A3) 
(A4) 
(A5) 



Note that Eqs. (|A3|) . (|A4|) and (|A5|1 still contain terms 
of 0{1/L), implicit in the shifts of the Fermi boundaries 
due to the creation of the high-energy particle as well 
as the low-energy charge and current excitations. Let us 
denote hy B — A(/+ = — /_ = A^/2) the Fermi boundary 
for the ground state of fermions. Then, 



5B+ = B+- B ^ 0{1/L), 
SB^ = B^+B ^ 0{l/L). 



(A6) 
(A7) 



To zeroth order in 1/ L, Eq. 
tion Eq. (lOi)) . 



3]) becomes the Lieb equa- 



We introduce the resolvent operator L formally by 

1 + • (l-^J (A,A') = 1, (A8) 



X(A|A')- 



2ti 



+B 



dy L{x\x")K{x" -y) 



K{X-X') 



2ti 



(A9) 



To zeroth order in Eq. (^3]) becomes Eq. (|6.3ip . 
with Pirap{X,X') = p.„-ap{X,X'\B,-B). By comparison 
with Eq. (jA9p . we have 



p,mp{X,X')=L{X\X'). 



(AlO) 
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The ground state energy (without the impurity) in the 
thermodynamic hmit is 



Egs — L 



+B 



dX eo (A) PCS (A) . 



(All) 



Consider a state with a single high-energy particle Xp 
and general Fermi boundaries B±. (The derivation is 
similar for the case of a single hole.) The energy is 



E = J2^o{Xj)+eoiXp), 



(A12) 



where the sum runs over all quantum numbers Ij be- 
tween /_ and /_(- . We expand Eq. (|A12p using the Euler- 
Maclaurin formula and obtain 



E 



L 



1 

24l 



dX! eo (A')p(A') + eo(Ap) 



p{B+) p{B-) 



(A13) 



The next step is to expand Eq. (jA13[) up to 0(1/L) 
using the expansion of p{\) in Eq. (jA2[) and assuming 
5B± - 0(1 /L). We find 



E 



E, 



GS 



e(Ap) 



TTV 

6L 



—plsiB) {L5B+Y + {L5B^ 



(A14) 



The third term on the rhs of Eq. (|A14[) is the stan- 
dard finite size correction to the ground state energy of 
a conformal field theory with central charge c — \. The 
difference between the energy of the excited state and the 
ground state energy is 

A£; = e {Xp) + ^plsiB) [{L5B+f + {LSB.f 

(A15) 

Now we evaluate the Fermi boundary shifts 6B± . From 
Eq. (|6.14p . we can write 



N 

T 



AN 



L 



dXp{X). (A16) 



with n^rnp and df^^ defined in Eqs. (jOO)) and (lOl) . 
The function Z{X) is the dressed charge defined by 



K 
2^ 



Z(A) = 1. 



(A20) 



The value of Z{X) at the Fermi boundary is related to 
the Luttinger parameter hy K = Z'^{B)}'^ The function 
^ (A) is defined by 



2e(A) = l- 



dA'L(A|A')- 



dA'L(A|A'). (A21) 



One can verify that the value of ^(A) at the Fermi bound- 
ary satisfies 



2Z{B)^{B) = 1. 



(A22) 



Using the definitions of Z{X) and ^(A), together with Eq. 
(jA10|) . we can write 



1 — 
imp 



„ = ±[Z(Ap,,0 - 1] 



(A23) 

2dP;^p = ±[2C(Ap,0-l]- (A24) 



Finally, substituting Eqs. IjAlSp and (|A19p into Eq. 
(jAlSp . we arrive at the finite size spectrum 



AE 



2ttv 



AN - rJ' 

"-imp 
2Z{B) 



-^Z^(B) {D-dlj\ 



(A25) 



APPENDIX B: PHASE SHIFTS IN TERMS OF 
SHIFT FUNCTION 

Here we prove the equivalence between our results^ 
in Eq. (|6.29p and (|6.30p and those of Cheianov and 
Pustilnik^^ in terms of the shift function. The shift func- 
tion F(A|A') is defined as the solution of the integral 
equationi^ 



Similarly, from Eq. (|6.15p . 



2D 
T 



dXp{X) 




1 
L 



El 



i j>/4 



dXpiX) 



(A17) 



Expanding the rhs of Eqs. (|A16p and (|A17p to 0(1/-L), 
we obtain 



AN 

IT 

D 

T 



L 

irnp 



{SB+ -5B^)pGs{B)Z{B), {A18) 
i6B+ + 6B_)pGs (S)^(B).(A19) 



Upon application of the integral operator on Eq. (IA2ip . 
we find that ^(A) obeys the equation 



from which we obtain the relations 



e(A) = ^+FiX\-B), 
e(-A) = ^^FiX\B). 



(B3) 
(B4) 
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Using Eq. (|A2ip . we can also show that 

e(A)+e(-A) = l. (B5) 
Hence from Eqs. (|A23P and (jA24jl . it follows that 



- F{X\-B)-F{X\B), 



(B6) 



2dLp = ^F{X\B)-F{X\-B). (B7) 



By taking the derivative of Eq. (jBl|) with respect to A 
and integrating either inside or outside the Fermi bound- 
aries we can show that 

ZiB)[F{B\X)^Fi-B\X)] = Z(A)-1, (B8) 
2aB){B)[F{B\X)+F{-B\X)] - 2e(A) - 1. (B9) 

We then obtain 



= -Z{B)[F{B\X)-F{-B\X)], (BIO) 



'imp 

dLp = -m[F{MB) + F{-B\X)]. (Bll) 



Finally, using Z{B) = {2£_{B))-^ = VK, we can write 
the phase shifts in Eq. (|6.33p in the form 



(B12) 



j^i^ = T27rVKFi±B\Xh). 



Eq. (jB12[) is equivalent to Eq. (14) of Ref. |16| for any 
finite value of B (away from half- filling). 



APPENDIX C: TWO-SPINON SPECTRAL 
WEIGHT FOR THE ZERO FIELD HEISENBERG 
MODEL NEAR g = tt 

In this appendix we determine what fraction of the 
total spectral weight is given by the two-spinon contri- 
bution to the dynamical structure factor in the zero field 
Heisenberg model of Eq. (|8.2p at g « tt and low fre- 
quencies. The calculation is based on a combination 
of Luttinger liquid and exact integrability techniques. 
Band curvature effects are completely ignored in the Lut- 
tinger liquid treatment. This appears to be valid near 
q = 2kF = for the case of zero magnetic field (half- 
filled band in the fermion model; see Sec. IVIIip . 

Based on Luttinger liquid methods it was argued^is^ 
that the asymptotic behavior of the equal time spin cor- 
relation function for the zero field Heisenberg model is 



G(a;,0) = (5|+,5|)(x 



i-irin'/'{\x\/a) 



(CI) 



Here a is a number of order 1. Later, the exact amplitude 
was determined^i^: 



Gix,0) 



{-ir\n'/\\x\/a) 
(27r)3/2|a;| 



(C2) 



The field theory methods indicate^''^ that G{x, t) is 
asymptotically a Lorentz scalar, depending only on x"^ — 



V t . This is true including the logarithmic corrections 
since they arise from a Lorentz invariant marginal term 
in the low-energy effective Hamiltonian. Therefore, the 
extension to the time-correlation function is 

i-lr\n'/V~v'{t-^ar]/a'] 
^ ' ' 47r3/2[x2-t;2(t-ia)2]i/2 ' ^^''> 

Here a > is of order 1/J. We now consider the Fourier 
transform, giving the structure function. Let us first of 
all ignore the log factor, setting it to one. Then we can 
do the integral exactly by changing variables to vt ± x 
and doing the two integrals separately, both of which are 
Gaussian. This gives: 



P (C4) 



S{q,uj) e{uj ~v\q~TT\) , 

-/tt^lj^ — v''(q — tt)^ 

for g w TT and small lu. We now include the log factor. 
S{q,uj) e{uj - v\q - tt\) . 

To see this we can change variables to 
(vt ± x~)(u}jv =F g) = u±, 

giving: 
S{q,uj) 



(C5) 
(C6) 



1 



87r3/2[t^2_„2(^_^)2]l/2 

gi(«++u_)/2 



du+du- 



X <^ In 



— («+ — ia)(u- — ia) 



a^[(Lu/vY-{q- Try 



1/2 



(C7) 



We now Taylor expand in powers of ln(— u+u-), integrat- 
ing term by term. The leading term gives Eq. (|C5p . 

To further justify this expression and understand an 
important subtlety, we compare S{t:+p), the equal time 
structure factor, at small obtained either by Fourier 
transforming Eq. (|C2p or by integrating Eq. (jCSP over 

UJ. 

Siq) = J dxe-'^-G(x) - J ^Siq,u). (C8) 
The first approach gives: 



Siir+p) = 



1 

(27r)3/2 
2 

(27r)3/2 



dx- 



\x\>a 



dxcosipx) ^ ' ' .{C9) 



Note that this integral diverges at |a;| — > so it was nec- 
essary to introduce a cutoff on the integration range. We 
have chosen it to be a for convenience but don't expect 



29 



the leading behavior at p ^ to depend on this choice. 
We now approximate this integral by 



S{Tr+p) 



(27r)3/2 7, 



dx ^ ' (CIO) 



where c is a constant of order 1. This is based on the 
approximation that, for small p, cospx is approximately 
constant and equal to 1 out to a large value of x of order 
The integral can then be done exactly giving: 



3(27r)3/2 



\n^/\c/{a\p\)] 



(Cll) 



(This result was obtained in Ref. l53l.) 

Now consider the second approach. Using uj^ ~ v\q ~ 
7r|, and introducing the high-energy cutoff, Z) cx 1/a cx J, 
we obtain 



S[n+p) 



2ti 



2 fVD'--'l duj' In^/^D/uj') 
dm' \n^l^{D/uj') 



duj' 



2ttuj' 



j\n'/\D/u;') 



3(27r)3/2 



(C12) 



ujI, in 



Here we have changed variables to uj' = 
the first step, extended slightly the upper cutoff in the 
second and used the fact that the denominator can be 
approximated by lo' down to a value of order at which 
we may simply cut off the integral, in the third step. 
The fact that Eqs. (|C11|) and (|C12|) agree gives further 
confidence in these calculations. 

The factors of 2/3 that arise from the integrals in both 
approaches are rather subtle and easily missed by more 
crude approximations. In evaluating the integral of Eq. 
(|C9p . for example, it is tempting to follows the same pro- 
cedure as we used in evaluating the integral in Eq. (jC7[) . 
changing integration variables Xo u = xp and approx- 
imating (Inu — lnap)^/2 ~ ln^/^(l/ap). However, this 
would be incorrect, giving an answer larger by a factor 
of 3/2. The fallacy with this approximation lies in the 
ultraviolet divergence of the integral at small u of order 
ap. Since the integral receives important contributions 
from this range of u it is not valid to drop the In u term 
inside the square root. The reason why we could make 
this approximation in evaluating Eq. (jC7p is that the 
integral is ultra-violet finite. 

Now let us compare Eq. (jCSP to the exact two-spinon 
result in Ref. After correcting typos, precisely the 
same asymptotic behavior, Eq. (jCSP is obtained, at g ~ 



TT and w — > except that it is smaller by a factor of 
\fC where \fC = 0.849829. Remarkably, the two-spinon 
approximation has the same combination of power-law 
and logarithmic singularities as does the exact result, but 
a smaller amplitude. 

Ref. 113 attempts to derive a formula for S{q) in Eq. 
(6.15). However, the tricky factor of 2/3 from the w in- 
tegral discussed above is missed. The correct two-spinon 
approximation to S{q) is 



(C13) 



smaller by a factor of vC w 0.85 than the exact answer. 
The fact that the two-spinon approximation to S((i) at 
q K, -K underestimates the correct answer by exactly the 
same factor, 0.849829, as does the two-spinon approxi- 
mation to S'(g, Lo) for q k, ix and small lo is not surpris- 
ing. It merely indicates that the dominant contribution 
to the frequency i nteg ral is from low frequencies. As was 
discussed in Ref. |2J, the total two-spinon intensity, in- 
tegrated over all q as well as all uj underestimates the 
exact answer (1/4) by a factor of 0.7289. So, we see that 
the two-spinon approximation does somewhat better at 
q near tt than at other values of q. 

Finally, we remark that the asymptotic behavior of the 
finite size equal time correlation function, with periodic 
boundary conditions, can be written^^: 

(27r)-'/^(_L/7r) sm(7r|a;|/L) 
The q = TT finite size equal time structure function: 



SiTT,L) 



, \n^/^\(L/TTa)sm(TTx/L)] 
dx- 



(L/tt) sin{TTx/L) 



(C15) 

can be approximated as in Eq. (jClOp . using the small x 
approximation to the integrand with some large x cut-off 
of 0(L): 



S{7T,L) 



dx 



(2^)3/2 X 



\n^/^{x/a) 



(27r)3/2 



ln='/2(L/a). 



(C16) 



Approximately this result was obtained^ by DMRG cal- 
culations, prior to the exact results, with l/(27r)3/2 ^ 
0.0632936 . . . replaced by 0.06789. We may obtain the 
two-spinon estimate of this quantity by replacing |g — 7r| 
by c/L in Eq. (|C13p . Again this is smaller than the 
exact result by the same factor of 0.85 and thus smaller 
than the old DMRG result by 0.91 rather than being dis- 
turbingly larger than the DMRG result as was stated in 
Ref. l2J, apparently due to missing the 2/3 factor. 
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